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Abstract 

A consistent analysis of relativistic pionic correlations and meson-exchange currents for electroweak 
quasielastic electron scattering from nuclei is carried out. FuUy-relativistic one-pion-exchange electromag- 
netic operators are developed for use in one-particle emission electronuclear reactions within the context of 
the relativistic Fermi gas model. Then the exchange and pionic correlation currents are set up fully respect- 
ing the gauge invariance of the theory. Emphasis is placed on the self-energy current which, being infinite, 
needs to be renormalized. This is achieved starting in the Hartree-Fock framework and then expanding 
the Hartree-Fock current to first order in the square of the pion coupling constant to obtain a truly, gauge 
invariant, one-pion-exchange current. The model is applied to the calculation of the parity-conserving (PC) 
and parity-violating (PV) inclusive responses of nuclei. Interestingly, in the pionic correlations terms exist 
which arise uniquely from relativity, although their impact on the responses is found to be modest. 

PACS: 25.30.Rw, 14.20.Gk, 24.10.Jv, 24.30.Gd, 13.40.Gp 

Keywords: Nuclear reactions. Inclusive electron scattering, Pionic correlations. Meson-exchange currents, 
Relativistic Fermi Gas. 



MIT/CTP#3167 

^Corresponding author. E-mail address: amaroQugr.es, phone: 34-958-246172, FAX: 34-958-249487 



Contents 



1 Introduction 



2 Parity-conserving electron scattering 

2.1 General formalism 



2.2 Pion exchange and correlation currents 

2^.1 Feynman diagrams and two-body currents) 
2.2.2 Particle- hole matrix elements 



2.3 Hartree-Fock renormalization in nuclear matter 



2.3.1 Nucleon propagator 

2.3.2 Nucleon spinors 

2.4 Self-energy current to first order 

2.5 Gauge invariance of the theory 

2.6 Nuclear hadronic tensor and electromagnetic response functions 

2.6.1 MEG and vertex pionic contributions! 

2.6.2 Relativistic self-energy responses 

2.7 Analysis of results 

2.7.1 MFC 

2.7.2 Correlations! 



3 Parity-violating electron scattering 

3.1 General formalism 



3.2 PV response functions 

3.3 Results 



4 Non-relativistic reductions 



4.1 Pion exchange currents! 

4.1.1 The r]F expansion 

4.1.2 Results for the responses using the relativized MFG model 

4.1.3 Comparison with the traditional relativistic corrections . . 

4.2 Pionic correlations 



5 Conclusions 



A Gauge invariance of two-body currents 



B Gauge invariance of the two-body current p-h matrix elements 



C Polarization propagator with nucleon self-energy 



D Renormalized self-energy response using the polarization propagator 



E On-shell self-energy and field strength renormalization function 



F The electromagnetic current operator 



89 
90 
93 
97 
100 
102 



2 



1 Introduction 



In modern experimental studies of electron scattering from nuclei [Q, ^] the typical values of 
energy and momentum transfer are comparable to or even larger than the mass scale set by the 
nucleon mass and accordingly one must expect relativistic effects to be important. Unfortunately 
the wave functions and operators used to describe this high-energy regime have been treated only 
approximately. Indeed it is still the case that many calculations continue to be done at the non- 
relativistic level with leading-order expansions of the electroweak currents involved |Q, |5|, ^, . 
However, a number of studies in recent years show that such an approach is highly constrained to 
work only at relatively low energies and momenta. 

In order to gain insight into which ingredients can or cannot be non-relativistically approxi- 
mated we have employed a simple model in which Lorentz covariance and gauge invariance can 
be maintained, namely, the relativistic Fermi gas model (RFG). Since our focus is placed on the 
quasielastic region where high-energy knockout of nucleons is kinematically favored, we believe 
that this model, while undoubtedly too simple to encompass the aspects of nuclear dynamics is 
nevertheless a convenient place to start in such explorations. Indeed, the problem of relativity 
in electroweak studies of nuclei is so difficult |9|, [l^, |Tl|, |2|, |l5|, |l8| that only in 

special frameworks such as the RFG can we hope to carry out all but rather severely approximated 
modeling. 

In the quasielastic regime we also expect pions to play a role that differs from the dynamics 
typically occurring near the Fermi surface where one expects other mesons (cr and uj in particular) 
to dominate. For quasielastic scattering the residual interaction of relevance is principally that 
between a low-energy hole and a very high-energy particle, and for this the pion is expected to 
play an important role. Accordingly, as the next step after the basic relativistic Fermi gas of non- 
interacting nucleons we have concentrated on one-pion exchange (OPE) effects in our description 
of the nuclear responses. These occur as correlation effects and also as two-body meson-exchange 
current effects. After developing approximate methods for modeling to this order, in recent work 
we have reached the stage where large classes of effects can be incorporated fully-relativistically. 
The present paper is a comprehensive discussion of what we have learned to this point based on 
this type of approach, together with comments on what directions future studies could follow. 

In particular, in a recent paper [^] we investigated the role played by pions in inclusive electron 
scattering from nuclei within the context of one-particle one-hole (Ip-lh) excitations, i.e., for the 
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dominant modes in the quasielastic regime. There we extended our previous work [20, 21| where a 
systematic investigation of relativistic effects in the nuclear electromagnetic responses spanning a 
wide range of kinematical conditions and accounting for both meson-exchange and isobar currents 
was carried out. In these studies a consistent first-order operator, embodying all Feynman diagrams 
built out of nucleons and pions with one exchanged pion and one photon attached to all the possible 
lines was set up to represent the two-body current. Importantly, the latter has been explicitly proven 



to be gauge invariant in |19|. 

In addition to the usual contact and pion-in-flight meson-exchange currents (MEC), this fully- 
relativistic operator includes as well the so-called correlation currents. The latter are often not 
included in model calculations because they give rise to contributions assumed already to be ac- 
counted for (at least in part) in the initial and final nuclear wave functions |^, However, 
our model is based on an uncorrelated relativistic Fermi gas whose states are Slater determinants 
built out of (Dirac) plane waves. Within a perturbative approach we are free to consider the one- 
pion correlation contributions to the responses as arising either explicitly in the wave functions 
or from an appropriate current operator acting on unperturbed states: our choice has been the 
latter. Clearly, should it be possible to sum up the whole perturbative expansion, then the results 
obtained starting with the true "correlated" wave function would be exactly recovered. 

In this paper we provide a deeper analysis of the impact of pions on the nuclear electromagnetic 



response in the Ip-lh channel |24, 25, 26, Just as for the MEC, the two-body correlation 

current also contributes in this sector and, of course, it should do so consistently, namely fulfilling 
gauge invariance at the level of one-pion-exchange. 

When the operator associated with a two-body current acts on the RFC ground state in general 
it changes the quantum numbers of two nucleons: the Ip-lh matrix element is then obtained via 
the integration of a one particle state over the Fermi sea. In the case of the correlation current two 



contributions are thus obtained. The first one, sometimes referred to as a vertex correction |29], 
arises from the exchange of a pion between the particle and hole; the second relates to the Fock 
self-energy (SE) |29, 3C] and dresses the particle and hole propagation lines. This one diverges. 



since it corresponds to a SE insertion on an external line, which field theory |31, |3^ tells us not 
to include in a perturbative expansion. Instead one should apply a renormalization procedure to 
dress the external lines by summing up the entire perturbative series of self-energy insertions. In 
the nuclear context this procedure leads to the relativistic Hartree-Fock (HE) approach. 



4 



In some relativistic calculations [29, 30| this contribution has been treated by introducing from 
the outset a Hartree-Fock propagator in the medium, which accounts for the SE diagrams. In- 
medium form factors for the Ip-lh current were also introduced neglecting however any momentum 
dependence in the self-energy and effective mass. Thus in [^0| the self-consistent Hartree mean field 
was inserted into the single-particle propagator, automatically including the Pauli blocking of NN 
pairs, whose contribution was thus included in the random-phase-approximation (RPA) responses 
computed there. A similar semi-phenomenological treatment of the nucleon self-energy in the 
medium at the non-relativistic level was implemented in [0]. More recently, the Dirac structure 



of nucleon self-energy in nuclear matter has been studied in |33], while a finite nuclei calculation 
based on the a-co model can be found in Q, where the relativistic Hartree model of |]l^] is used for 
the single particle bound states. 

In the difficulty of the SE insertion in first order was avoided by computing the associated 
self-energy response as the imaginary part of the corresponding polarization propagator with one- 
pion-exchange SE insertions on the particle and hole lines. A finite result was thus obtained in first 
order (one pionic line) without resorting to the HE approach. The question then arises whether 
it is possible to obtain the same result for the self-energy response function starting with finite 
well-defined matrix elements of the current operator. 

In this paper we answer this question by constructing a renormalized self-energy current cor- 
responding to one-pion-exchange. This current acts over free Dirac spinors and leads to the same 
response functions as those obtained by taking the imaginary part of the polarization propagator 
computed to first order. It should be clear that in this work the concept of renormalization has a 
many-body significance, namely it amounts to a relativistic HE approximation and ignores (see ||l9[| ) 
the additional vacuum renormalization due to the change of the negative-energy sea induced by 
the nuclear medium |]30t| . 

The new current is obtained by renormalizing spinors and energies and by expanding the result- 
ing in- medium one-body current to first order in the square of the pion- nucleon coupling constant, 
to be consistent with the requirement of dealing with diagrams having only one pionic line. The 
renormalized quantities should then be obtained in the general case by solving a set of self-consistent 
relativistic HE equations numerically. However, one of the goals of this paper is to show that to first 
order the solutions and the corresponding corrections to the bare single-nucleon current operator 
can be expressed analytically in terms of a simple electromagnetic operator. This operator accounts 
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for two main effects induced by the interaction of the nucleon with the medium: the first is the 
enhancement of the lower-components of the Dirac spinors; the second is a global renormalization of 
the spinors in the nuclear medium. These effects are genuine relativistic corrections that are absent 



in a non-relativistic framework [24|. Actually a third renormalization effect also arises, related to 
the in-medium modification of the energy- momentum relation for a nucleon, which is here treated 
in first order of the square of the pion-nucleon coupling constant (in other approaches this effect is 
embedded in a constant effective mass p9[). 

Using the renormalized SE current operator together with the MEC and the vertex exchange 
operator we prove the full gauge invariance of the current if account is taken of the change in 
energy arising from the HE renormalization to first order. The results for the inclusive response 
functions we obtain with this current agree completely with the ones of jl^, where the polarization 
propagator technique was used. 

The present review is organized as follows: in Section 2 we focus on parity-conserving electron 
scattering from nuclei. We begin in Section 2.1 with some general formalism and then in Section 2.2 
discuss the pion exchange and correlation currents. There we revisit the full set of Ip-lh current 
operators with one pion-exchange line which contribute to the electro-excitation process, paying 
special attention to the SE contribution. We show the necessity of re-defining the otherwise infinite 
self-energy diagrams. In Section 2.3 we develop the Hartree-Fock renormalization scheme as a 
vehicle to addressing this problem, going on in Section 2.4 to expand the renormalized spinors and 
energies to first order in the pion coupling constant squared obtaining a new self-energy current. 
Then in Section 2.5 we prove the gauge invariance of the theory. To conclude this section we go 
on to discuss the hadronic tensor and electromagnetic response functions (Section 2.6) and present 
some typical results (Section 2.7). In Section 3 we briefly discuss parity- violating electron scattering 
to place it in context with the above studies. In Section 4 we make contact with non-relativistic 
expansions schemes, both for the pion exchange currents (Section 4.1) and for the pionic correlations 
(Section 4.2). In Section 5 we summarize our results and draw our conclusions and end with a series 
of Appendices where more technical aspects of the formalism are compiled. 



2 Parity-conserving electron scattering 
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2.1 General formalism 

The general formalism involved in the description of (e,e') processes for quasielastic kinematics has 



been derived and discussed at length in several papers (see for instance |3^, ^ |3 



40, Here we summarize only those aspects that are of special relevance to the analysis that 

follows. We limit our attention to the Plane Wave Born Approximation (PWBA), i.e., the electron 
is described as a plane wave and interacts with the nuclear target via the exchange of a virtual 
photon. The laboratory system variables involved in the process are = (e, k) and K''^ = (e' , k'), 
the initial and scattered electron four-momenta, and P-^ = (£'j,pj) = (Mj,0) and Pj = {Ef,pf), 
the initial and final hadronic four-momenta, respectively. The four-momentum transferred by the 
virtual photon is Q'^ = {K — K')^ = {Pf — Pi)^^ = (w,q); for electron scattering the momentum 
transfer is spacelike, Q"^ = uj"^ — < 0, with q = |q|. The 5-matrix element in PWBA can then be 
written as 

5/, = -27ri5(£;/ -S, -a;)^(k',s'|ieM(0)|k,s)(/|J>(g)|i), (1) 

where 

(k',/|je;.(0)|k,.) = (^^^y^'^,,(k')7^n.(k) (2) 

is the electron current matrix element and J''((5) is the Fourier transform of the nuclear electro- 
magnetic current operator. 

We assume Lorentz invariance, parity conservation and work in the extreme relativistic limit 
(ERL), in which the electron energy e ^ mg. Under these conditions the unpolarized, inclusive 
(e,e') cross section reads 
da _ 2o? fe'^ 

Here a is the fine structure constant and Jig the scattered electron solid angle. The term o"m 
represents the Mott cross section which in the ERL reduces to 

2 



^rji.uW'"' = aM[vLR^iq,uj)+VTR^iq,oj)] . (3) 



a cos 6e/2 
2e sin^ ee/2 



where Bg is the electron scattering angle, and rj^u and W^'^ are the leptonic and hadronic tensor, 
respectively. Within PWBA the leptonic tensor simply reads 

^^, = K^Kl + K'^K, + ^g^,. (5) 
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The kinematic factors vl and vt are evaluated from the leptonic tensor using standard techniques 
(see, for example, 



whereas the longitudinal and transverse (with respect to the momentum transfer q) response func- 
tions and are constructed directly as components of the hadronic tensor W'^'^ according 
to 



i?^(g,u;) 



q 



(8) 



R^{q,Lo) = + (9) 

where we use a coordinate system with the z-axis in the direction of the vector q. Note that if 
gauge invariance is fulfilled, implying that W'^^ = W^'^ = {uj/q)W'^^ and W^^ = {uj / q)'^W'^^ , then 
R^ is simply the time component of the hadronic tensor, namely W^^. Hence R^ is determined by 
the charge distribution, whereas R^ reflects the current distribution of the nuclear target. 

The hadronic tensor and consequently the response functions derived from it embody the entire 
dependence on the nuclear structure, specifically on the charge and current distributions in nuclei, 
and accordingly these provide the prime focus in analyses of electron scattering. There are various 
options on how to proceed in performing such analyses (see, for example, ^^), depending on the 
specific problem under consideration and on the approximations to be made. In what follows 
we recall two common expressions for the hadronic tensor W^"^ and comment briefly on their 
applications. 

First, the hadronic tensor can be defined according to 

i f 

where J^{Q) represents the nuclear many-body current operator, the nuclear states \i) and |/) are 
exact eigenstates of the nuclear Hamiltonian with definite four-momenta, and the sum with a bar 
means average over initial states. This form is very general and includes all possible final states 
that can be reached through the action of the current operator J^{Q) on the exact ground state. 
In our perturbative approach we shall use eigenstates of the free Hamiltonian Hq (which describes 
the free relativistic Fermi gas) and include correlations among nucleons in the current mediated 



by the exchange of pions. This current of course allows one to reach both the p-h and the 2p-2h 
sectors in the Hilbert space of Hq. In the present work, however, we shall restrict our attention to 
the former. 

A different option for evaluating the nuclear responses exploits the polarization propagator 11^'^ 
(also referred to as the current-current correlation function). The latter can be expressed in terms 
of the full propagator, G, of the nuclear many-body system, since closure can be used to carry out 



the sum over the final states in eq. (10). Then one has for the hadronic tensor |43] 



p^M- = _ilnan'^-(g,g;c^) =-ilmV(i|jtM(Q)G(^ + S.)j-(Q)|i) . (11) 
TT vr ^—^ 

A possible advantage of this approach relates to the existence of a well-defined set of rules (the 



relativistic Feynman diagrams) which allows one to compute H^'^ perturbatively |43]. 

Obviously the two procedures are equivalent and hence the observables calculated using the 
expressions for the hadronic tensor given by eqs. (|lO|) or (|ll|) should be the same. However, notice 
that eq. (|l^) is less suitable for dealing with situations where the nuclear current matrix element 
{f\J^\i) is divergent. In this case one proceeds either by computing directly the responses via the 
polarization propagator or by first renormalizing the current matrix element and then by using 
eq. (P. 

Finally, we remark that gauge invariance must be fulfilled both at the level of the nuclear 
current matrix elements and at the level of the hadronic tensor and/or the polarization propagator. 
A consequence is that the electromagnetic continuity equation should be satisfied. In other words 
in momentum space all of the expressions Q^{f\J^{Q)\i), QfjW^^ and Qfjl^^ should vanish. 

2.2 Pion exchange and correlation currents 

Working within the framework of the relativistic Fermi gas (RFG) model, i.e., for nucleons moving 
freely inside the nucleus with relativistic kinematics, in this section we present a detailed study of the 
electromagnetic currents accounting for the effects introduced by pions in first-order perturbation 
theory (one-pion exchange). 

2.2.1 Feynman diagrams and two-body currents 

The linked, two-body Feynman diagrams that contribute to electron scattering with one pion- 
exchange are shown in Fig. 1. The first three correspond to the usual meson-exchange currents 
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Figure 1: Feynman diagrams contributing to the two-body current with one pion-exchange. The 
wide line in the correlation diagrams (d)-(g) means a fully-relativistic Dirac propagator for the 
nucleon. 

(MEC): diagrams (a), (b) refer to the contact or seagull current, diagram (c) to the pion- in-flight 
current. The four diagrams (d)-(g) represent the so-called correlation current and are usually not 
treated as genuine MEC, but as correlation corrections to the nuclear wave function. However, 
again we note that our approach puts all correlation effects in the current operator and uses an 
uncorrelated wave function for the initial and final nuclear states. 

In this work we use Bjorken and Drell conventions and pseudovector ttNN coupling (the 
effect of a pseudoscalar coupling will be commented on later), namely 



/ _ 



(12) 



where ij} is the nucleon field, (pa is the isovector pion field, / represents the ttNN coupling constant 
and rriTT is the pion mass. The electromagnetic currents corresponding to diagrams (a)-(g) are 
obtained by computing the ^-matrix element 



^fi — Sfi[P[.,P2-,Pl-,P2) - Sfi{P[,P2,P2,Pl] 



(13) 
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for the absorption of a virtual photon by a system of two nucleons, namely for the process 

^ + Ni+N2^ N[ + N^ , (14) 

with Pi, P2 {P{, P2) being the initial (final) four-momenta of the two nucleons involved (see Fig. 1). 
The electromagnetic current is then defined according to 

SfiiPi,Pi,Pi,P2) = -ieAi,{Q)27r6iE[ +E'^-E^-E2- uj){P[P^\f{Q)\PiP2), (15) 

where Af^{Q) is related to the matrix element of the electromagnetic field between the incident 
photon with momentum Q and the vacuum state, namely 

(0|^^(X)|7(Q)) = ^^(g)e-«-^. (16) 

Finally, the on-shell matrix element of the two-body current can be written in terms of a function 
(p'l , P2 > Pi > P2 ) as follows 

{pip^\fmpiP2) = 

= (27r)^<53(p;+p-,-q-pi-p2) ,,,, j'{p[,p'2,Pi,P2), (17) 

V'^[Ep,Ep2Ep>^Ep0'-/'' 

where m is the nucleon mass, V is the volume enclosing the system and Ep = i/p^ + the 
on-shell energy of a nucleon with momentum p. The four-momenta — indicated by capital letters 
— are defined in Fig. 1. 

The general relativistic expressions for the seagull (diagrams a,b), pion- in-flight (c) and corre- 
lation (d-g) current matrix elements are (isospin summations are understood) 

• Seagull or contact: 

/2 _ pV 
(P'1>P2>P1>P2) = —ie3abuip'i)Ta'y5 ^in(pi)— 3 2'^iP2)Tb75l^u{p2) + (1^2) (18) 

• Pion-in-fiight: 

s2 jp (K\ 1^2)^ 

ip(pi,P2,Pl,P2) = —2ie3ab , 2 ^Imk"^^ ^(pOraTS ;^m(Pl)^(P2)7-675 ;^2u(P2) 

(19) 

• Correlation: 
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/2 _ 

Jcor(p'l'P2iPl'P2) = — 2^(p'l)Ta75 ;^l^^(Pl 



X ii(p^) [r„75 /^i5f(P2 + g)r^(g) + r'^(Q)5F(p^ - g)Ta75 ^i] u(p2) 

+ (1^2). (20) 

In the above, Ki, K2 are the four-momenta given to the nucleons 1, 2 by the exchanged pion, and 
they are defined in Fig. 1 , while and i<V are the electromagnetic isovector nucleon and pion form 
factors, respectively. Furthermore, Sf{P) is the nucleon propagator and ^^{Q) the electromagnetic 
nucleon vertex, i.e., 

Sf{P) = (21) 
r'^(Q) = FiT^^ + -^F2a'^-Q, , (22) 



Fi and F2 being the Dirac and Pauli form factors: for these we use the Galster parameterization |44]. 
Finally, the spinors (for brevity we denote n(p, Sp) by m(p)) are normalized according to the Bjorken 
and Drell convention and the isospinors are not explicitly indicated. 

The seagull and pion-in-flight currents shown above coincide with the expressions given by 
Van Orden and Donnelly |45| if account is taken for the different conventions used for the gamma 



matrix 75 and for the metric. Concerning the correlation current note that, at variance with |45], 
it embodies both the positive and negative energy components of the nucleon propagator. 

A crucial point to be stressed is that the sum of the relativistic seagull, pion-in-flight and 
correlation currents satisfy current conservation, i.e. Q^J^ = 0, provided some assumptions are 
made for the form factors involved in the various currents. This is proven in Appendix A (see 
also |l^) where we show that when the seagull and pion-in-fiight currents are multiplied by the 
same electromagnetic form factor Fi , gauge invariance is fulfilled, i.e. 

Qf^ij^+jiS+jD = , (23) 



where the two body currents in eq. (p^) are defined in eqs. (18-20). 

It is also possible to use different phenomenological electromagnetic form factors for the 
nucleon and pion — even introducing phenomenological form factors at the strong pion-nucleon 
vertices — without violating current conservation, by appropriate modification in the currents 
through the generalized Ward-Takahashi identity Wi 
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2.2.2 Particle-hole matrix elements 

In this report we deal with the case of one-particle emission induced by the two-body currents 
introduced above. The matrix element of a two-body operator between the Fermi gas ground state 
and a Ip-lh excited state reads 

= E ' (24) 

k<F 

where the summation runs over all occupied levels in the ground state, and thus includes a sum 
over spin (s^) and isospin (t^) and an integral over the momentum k. 

The first and second terms in eq. (p^) represent the direct and exchange contribution to the 
matrix element, respectively. It can be easily verified (see, e.g., |19, pH]) that in spin-isospin 



saturated systems the direct term vanishes for the MEC and pionic correlation currents upon 
summation over the occupied states. Hence only the exchange term contributes to the p-h matrix 
elements. The associated many-body Feynman diagrams are displayed in Fig. 2. Diagrams (a,b) and 
(c) correspond to the seagull (or contact) and pion-in-flight contributions, respectively. Diagrams 
(d-g) represent instead the correlation contributions. Here we distinguish the exchange of a pion 
between a particle and a hole line (d,e), giving rise to the so-called vertex correlation (VC), and 
the self-energy insertions on the nucleonic lines (f,g). After carrying out explicitly the sums over 
the internal spin, s^, and isospin, t^, the fully-relativistic expressions for the MEC (seagull and 
pion-in-flight) and correlation (vertex correlations and self-energy) currents turn out to be 

• Seagull 

• Pion-in-flight 



Jp^(P,h) 

''^'V^/ ' E^[(P-Kr- ml] [{K -HY-ml 



2m— legafe > T^TTT^ TTTo nTTTT? TTv) 7T^u{p)Tan[fi - m)u[h.) 



(26) 



• Vertex correlations 
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+ {P- - m)-i^raV^{Q)TaSF{K - Q)75 | n(h) 

= ^M + ^z^ (27) 



• Self-energy 



= W^+H;: (28) 



The effects of the medium are included through the summation in eqs. (25,28) over the intermediate 
momentum k up to the Fermi momentum. 

In the thermodynamic limit the sum y '}2]i<kp becomes an integral over the momentum / 
in the range < A; < kp^ kp being the Fermi momentum, and over the angular variables 6k, (pk- 
Note that, although the global factor ——=== has been extracted from the current in eq. (|2^ , 

TTl 

the factor , associated with the internal line, has to be retained inside the sum. Note also 
that, in order to fulfill gauge invariance, we have assumed F.,^ = FY ■ 

The vertex p-h matrix element splits into two terms and representing the forward- and 
backward-going contributions, respectively (Fig. 2d and 2e). They are 

= -V-2 E ^HPHiP- P)Sf{K + Q)r,r>^{Q)raj,^^^—^±^u{h) (29) 
V^i k<fcV [K - HY - mi 

B' = -T^i: ^^iP) (p ^J^r 2 75Tar^(Q)ra5i.(j^ - Q)l,{P- ^)n(h) . (30) 

Similarly, the self-energy p-h matrix element splits into two terms, and Tiy^. The former 
corresponds to the diagram with the pion inserted in the particle line (Fig. 2g), whereas the latter 
describes the diagram with the pion inserted in the hole line (Fig. 2f). They are given by 

K = f'^ipx A- - m)j^I0-^^s,(P)m,QH^ (31) 

K = 1^ g^lPjr-lQj^Hg) %!'l„^ (f^ - . (32) 
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Figure 2: Many-body Feynman diagrams contributing to the one-body current with one pion- 
exchange. The thick Une in the correlation diagrams (d)-(g) corresponds to a fuUy-relativistic 
Dirac propagator for the nucleon. Diagrams (d)-(e) represent the vertex current, while diagrams 
(f) and (g) represent the self-energy current of the hole and of the particle, respectively. 



Finally, splitting also the electromagnetic nucleon operator into its isoscalar and isovector 
parts, one obtains the isoscalar and isovector contributions to the self-energy and vertex p-h matrix 
elements. The final results can be cast in the form 

= -^^E £^(P)7.(A-- P)S.(K + Q)rM^>(Q),. u(h) (33) 

for the isoscalar and 



X Sf{K + Q)r^(^)(Q)(r3 + iesabTanH 

f2 



{f[ - m) 



ru(h) 



{K - Hf - ml 

Vmi El, [P - Ky - mi 



(35) 



(36) 
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for the isovector vertex p-h matrix elements, and 



up 



(37) 



and 



(38) 



for the self-energy matrix elements. Interestingly, the isoscalar /isovector ratio is = —3 in the 
vertex matrix element, whereas in the self-energy case it is the unityQ. Note that the MEC (pion- 
in-flight and seagull) p-h matrix elements are purely isovector, whereas the vertex and self-energy 
correlations get both isoscalar and isovector contributions. 

The VC and SE p-h matrix elements involve the nucleon propagator Sf{P) which in some 
situations may imply the occurrence of singularities. In the case of the vertex diagrams, the four- 
momenta appearing in the propagators are K + Q and K — Q for the forward- (Fig. 2d) and 
backward-going (Fig. 2e) contributions, respectively, and an integration over k should be done. For 
q > 2kp (no Pauli blocking) it can be proven (see ||l^) that only the forward diagram contains a 
pole, i.e., a value of the inner momentum k exists such that the nucleon carrying a four-momentum 
K + Q is on-shell. In this situation the forward vertex p-h matrix element is evaluated by taking the 
principal value in the integral over cos 0^ . In the case of the backward-going diagram the nucleon 
propagator Sp{K — Q) has no singularity for the kinematics in which we are interested. 

The case of the self-energy diagrams is clearly different. Here the particle {p) and hole (/i) 

are described in the Fermi gas by unperturbed plane waves, i.e., they are on-shell, and hence the 

propagators Sf{P) and Sf{H) diverge. The divergence of the diagrams (f)-(g) is reminiscent of 

the well-known infinity occurring in standard perturbative quantum field theory, when self-energy 

insertions in the external legs are included in Feynman diagrams [32|. As is well-known, there one 
^ The latter result stems from the relation ts + ie3abTaTb ~ — ts; however we prefer to leave the isospin structure of 
the isovector exchange as in eqs. (^sj^), since it makes more transparent the self-energy and exchange cancellation 
in the continuity equation, as shown in Appendix B. 
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Figure 3: Diagrammatic representation of the self-energy current for a hole (a) and a particle (b). 
The self-energy is defined to first order as the Fock insertion shown in (c) with one pion-exchange. 

should renormalize the theory by dressing the external legs, propagators and vertices. In the nuclear 
matter case we assume that the particle-physics effects are already accounted for by the physical 
masses and electromagnetic form factors. However, an additional nuclear physics renormalization, 
arising from the interaction of a nucleon with the nuclear medium, should be included at the 
one-pion-exchange level to account for the self-energy diagram. 

The self-energy current in eq. (|28|) can be written in the following form: 

j^^(p,h) = uip)E{P)SF{P)T^'{Q)u{h)+uip)T^^iQ)SF{H)E{H)u{h) , (39) 

where is the nucleon self-energy matrix that in first order reads 

^ ]i<kp Sfc.tfe \ J TT 

This is diagrammatically displayed in Fig. 3. The SE matrix, shown in Fig. 3(c), corresponds to 
the Fock term of the mean-field potential (the Hartree or direct term is zero for pion exchange, 
since it involves a pion carrying zero momentum). 
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Performing the sum over the internal spin (s^) and isospin (t^) indices and using the com- 
mutation properties of the gamma matrices to ehminate 75, the self-energy can be cast in the 
form 

where the sum over k has been converted into an integral. Note that the self-energies and 
appearing in eq. ( |39|) are evaluated for free particles and holes, i.e., for P'^ and on-shell. 
Hence the self-energy contributions to the current are divergent, since so are the free propagators 
Sp{P) and Sp{H) in eq. (^). Therefore they should not be computed using eq. (p9[), but rather 
one should first renormalize the wave function and the propagator of the particles in the medium. 
This is achieved through the summation of the full series of diagrams with repeated self-energy 
insertions displayed in Fig. 4. 

Now the energy of a particle in nuclear matter is modified by the interaction with the medium 
and, as well, through its energy-momentum relation. Thus the associated momentum is no longer 
on-shell and therefore in the next section we shall evaluate the self-energy for off-shell particles. 
In the first iteration, corresponding to one pion-exchange, the particle is off-shell, but the 
intermediate interacting hole K^^ is still on-shell0. In this case, with the help of Dirac spinology, 
one writes 

{P- - m){P- fC) = 2{P ■ K - m^){P + m) - (P^ - m^){f<[ + m) , (42) 



which allows one to recast the self-energy in eq. (|4l| ) for the off-shell momentum P in the form 

^ ' ml J (27r)3 ^ ^ '2Ei, {P - Kf -ml ^ ' 

Note that the second term inside the integral vanishes for P on-shell. 

In general the self-energy of a nucleon in nuclear matter can be written in the form [pS| ]: 

= mA{P) + B{P)-foP^ - C{P)^ ■ p . (44) 

In contrast to the quantum field-theory decomposition S(P) = mA + B /P, owing to the non- 

invariance under a boost of the step function 0{kp — k) appearing in the self-energy, in nuclear 

matter B{P) 7^ C{P). This in turn reflects the existence of a privileged system, namely the lab 
^Note that in deriving eq. ( ^ ) we have assumed free spinors u{k)\ hence eq. (^) is only valid for K'^ on-shell. 
The off-shell case requires a redefinition of the spinors ?t(k) according to an interacting Dirac equation, as is shown 
later. 
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system where the Fermi gas has total momentum Pfg = 0- Here it is natural to compute the 

self-energy. Under a boost, the Fermi gas ground state is no longer characterized hy k < kp and 

also the self-energy takes a different form. 

In the case of the Fock self-energy in eq. ( ^3|) the functions A,B,C can be expressed in terms 

of the integrals (for on-shell) 

d^k 1 1 

-^6{kF - k)- 



m ^ 

L^(P) = 
Indeed one gets^ 
A{P) = 
B{P) = 
C{P) = ■ 



(27r) 
j^Oikp - k) 



2Ei, {P - K f - ml 

1 Ki" 
2Ek {P - KY - ml ■ 



L 

m\ 

-3^<!2 

mi 

m± 



3^ {2 [Pf.L^'iP) - m^I{P)] - (P2 _ m^)I{P)} 



P^L'^iP) - m^I{P)] - (P2 _ ,„2)^o(^) 



Po 



Pf.L^'iP) - m^I{P) - (P2 _ m^) 



LsiP) 



P 



(45) 
(46) 

(47) 
(48) 
(49) 



Note that A = B = C for P on-shell. In this case one simply has S(P) 



on— shell 



A{P){m+ P). 



2.3 Hartree-Fock renormalization in nuclear matter 

In this section we discuss the renormalization of the nucleon propagator and spinors associated 
with the pionic self-energy in a Hartree-Fock scheme. 

2.3.1 Nucleon propagator 

The Hartree-Fock (HF) nucleon propagator in the nuclear medium is the solution of Dyson's equa- 
tion 



Shf{P) = So{P) + So{PMP)SHFiP) , 
where S(P) is the HF proper self-energy and 



(50) 



(51) 



Z' + m + ze P + m — iepo 

is the free propagator in the medium. Equation ( pO[ ) results from summing up a series with an 

infinite number of self-energy insertions (see Fig. 4) for each of the two terms in (|5l[)P| , namely 

^ L is parallel to p since, choosing p along the z axis, the azimuthal integration in eq. (^6|) yields Lx ~ Ly — 0. 
''No interference term arises, since 6{k — kF)6{kF — k) =Q. 
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Figure 4: Diagrammatic series for the nucleon propagator in the medium. 

^ +^r—^iP)-Fr— + ^r—^iP)^r—^iP)^r—+--- = -Fi ^ -(52) 



P — m P — m P — m P — m P — m P — m P — m — S(P) 

Using the spin decomposition of the self-energy in eq. ( |4^ ) , we can write 

p-m- T.{P) = [1 - B{P)\ 70P0 - [1 - C{P)] 7 • P - [1 + A{P)] m . (53) 

Now the new four-momentum = f^[P), which is related to P^ as follows 

/o(^') = \ I Po (54) 

f(P) = p, (55) 

and the functions 

HP) = I ^ ^[p] rn (56) 
m = (57) 



(58) 



allow one to recast eq. (|52| ) in the form 

1 _ z{P) _ z{P) 



P-m- S(P) 7o/o(P) - 7 • P - fh{P) f{P) - m{P) ' 
For a nucleon with a fixed three-momentum p, the pole of the propagator in eq. (^) in the variable 
Pq defines the new energy of the nucleon in the medium. To find the latter we introduce 



E{P) = E{p, m{P)) = ^p2 + m(P)2 . (59) 

Then the propagator reads 

1 _ z{P) f{P)+m[P) ^g^^ 



p-m- S(P) /o(p) _ E{P) fo{P) + E{P) 
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and its pole po is found by solving the implicit equation 



MP)=E{P) , 

which, exploiting eq. ( 
1 - C(P) 

Po - 



can be recast as follows 



l-C(po,p) 



p2 + 



p2 + m(po,p)^ • 



(61) 



(62) 



l_ij(P)V- ^ ^ l-B{po,p) 
The solution of eq. ( p^ ) for fixed p defines the new dispersion relation po = e(p) for interact- 
ing nuclear matter. Once the above equation has been solved, the field strength renormalization 
constant 

ziP)_ 

PO=e(p) 



Z2(p) = Res 



(63) 



/o(P) - E{P) 

defined as the residue of the first factor on the right-hand side of eq. (|6^) at po = e(p), can be 
computed. Indeed using eq. (|57|), ^2(p) is obtained by expanding the denominator around the pole 
e(p), i.e., 



[1-C{P)] UP)-E{P) =Z2(p)-MP0-e(p)] + --- ; 



(64) 



hence 



^2(P)-^ = ^ 

dpo 
d 

dpo 



PO=<:(p) 



{[l-C(P)] [fo{P)-E{P)]} 
{[l-B{P)]po-[l-C{P)]E{P)} . 



PO=e(p) 

With the help of eq. (]59|) the derivative can be easily evaluated, the result being 



l-B 



dB m dA dC 

dpQ E opo E opo 



(65) 



(66) 



Jpo=e(p) 



2.3.2 Nucleon spinors 



The self-energy modifies not only the propagator and the energy-momentum relation of a nucleon, 
but, as well, the free Dirac spinors. In fact the spinors are now solutions of the Dirac equation in 



the nuclear medium |49], i.e., 

[P-m-^{P)]<P{p)=0 , 
which, again using the decomposition in eq. (|4^), can be recast as follows 

hoUP)--fP-m{P)]<P{p) = 0, 



(67) 



(68) 
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the functions fo{P) and m{P) being defined in eqs. (^) and (p6|), respectively. Equation (pSD has 
the same structure as the free Dirac equation; hence for the positive-energy eigenvalue one has 

/2(P)=p2+m2(P) , (69) 

which implicitly yields the energy pQ = e(p) of the nucleon in the nuclear medium. This result was 
already obtained as the pole of the nucleon propagator. Then the corresponding positive-energy 
spinors (s being the spin index) read 



^(p) + m(p) 
2m(p) 



1/2 



Xs 

a ■ p 



E{p) + m(p 

where the two functions of the three-momentum p 



-Xs 



^2(p)ws(p,m(p)), 



(70) 



m(p) = m(e(p),p) 



(71) 



and 



Eip) ^ Eieip),p) = + m(p)2 , (72) 
represent the nucleon effective mass and effective energy corresponding to po = e(p). The field 



strength renormalization constant, y/ (p) , of the new spinors, defined in eq. (pq), is required by 
renormalization theory, since the propagator in eq. ( ^8|) for po close to the pole e(p) reads, from 
eq. (H, 

1 Z2(p) /(p)+m(p) 1 m(p) 



(73) 



P-m-Y.{P) po-e(p) 2S(p) Po - e(p) ^(p) 

Once the new spinors have been computed, the self-energy can be evaluated by inserting (/>(k) 
instead of u(k) into eq. (|40|). Then the Dirac equation should be solved again with the new self- 
energy and so on. This self-consistent procedure leads to the relativistic Hartree-Fock model which 
has to be dealt with numerically. 

In this paper we do not attempt to solve the HF equations, since we are interested only in the 
OPE first iteration correction to the single- nucleon current. Although the latter cannot be derived 
by directly applying the Feynman rules, it can still be identified with the self-energy diagrams of 
Fig. 2 (f)-(g). Thus in the next section we shall compute the renormalized one-body current using 
the new spinors and energy-momentum relation and then expand it in powers of the square of the 
pion-coupling constant /^/m^. As we shall see, the unperturbed one-body current is thus recovered 
as the leading-order term whereas the first-order term is the searched for self-energy contribution. 
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+ 




+ 




(a) 



(b) 





(d) 



Figure 5: Diagrammatic series for the one-body electromagnetic current with dressed external 
lines. 

It is also important to remark that the use of the new 'renormalized' wave functions (j)s leads to 



a slightly modified global momentum distribution as shown in |50|. Note however that the number 
of particles is conserved without modifying the value of the Fermi momentum selected. 



2.4 Self-energy current to first order 

The particle-hole (p-h) current matrix element in the HF approximation reads 



(74) 



where the spinors 0(p), the first iteration solution of the Hartree-Fock equation, are given by 
eq. ([70|). Hence eq. ( |7^ represents the electromagnetic excitation of the p-h pair with dressed 
external lines corresponding to the sum of the diagrams shown in Fig. 5. In order to obtain 
a genuine one-pion-exchange expression we expand eq. ( [7^ ) in powers of the square of the pion 
coupling constant /^/m^ and single out the first-order term, i.e. the one linear in f"^ /m^. We shall 
still refer to the current thus obtained, representing the OPE contribution, as the "self-energy" 
current and, importantly, we shall show that it yields a finite contribution, free from the divergence 
problem of the current in eq. (^) . 

To proceed we start by deriving the HF energy e(p) to first order in /^/m^. For this purpose 
we note that the functions A[P), B{P) and C{P) defined in eqs. are of order 0(/^/m^). 
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Hence the following expansion of the Dirac mass in eq. (|5q ) holds: 

^^^'^ = "^ I-cIp] = + ^^^^ + + ^ (^) • ^^^^ 

Inserting this into eq. (|6^ ) for the energy and expanding again to first order in /^/m^, we get 



= e(p) ~ [1 - C{P) + S(P)] ^p2 + m2 + 2m2 [^(P) + C7(P)] 

= Ep + AE{po,p) , (76) 



where Ep = -y/p^ + rr? is the unperturbed free energy and 



1 



AE{po, P) ^ — m'A{P) + E^B{P) - p'C{P) +0[^] (77) 



/ 



4 



m. 



4 



is the first-order correction to the energy. Next we can insert the above value of pq inside the 
argument of the functions A, B, C. Expanding the latter around the on-shell value po = we get 

A{P) = A{po, p) = A{Ep + AE, p) = A{Ep, p) + O (j^^^ Ao{p) , (78) 

where ^o(p) = ^{Ep,p)- Likewise to first order we obtain 

B{P) ^ B{Ep,p) ^ Bo{p) (79) 
C{P) ^ C(Pp,p)^Co(p) . (80) 

Recalling that for P on-shell the functions A, B, C coincide, i.e., ^o(p) = Bq{p) = Co(p), we can 
insert these on-shell values into eq. (77) and, neglecting terms of second order, i.e., O we 
finally arrive at the result 

= e{p) = Ep + ^Bo{p){m' + El-p') + I ^] = Ep + ^2m' Bo{p) + I ^] .(81) 

The above expression can be recast in terms of the on-shell value of the self-energy 

So(p) = 2mSo(p) , (82) 
which satisfies the relation 

S(£;p,p)n(p) = Eo(p)u(p) , (83) 

thus showing that the free spinors are eigenvectors of the on-shell self-energy matrix S(£'p,p) 
corresponding to the eigenvalue So(p). Hence to first order in f"^ /rn^, the HF energy in eq. ( |8l] ) is 
found to read 

e(p):^Pp + ^So(p) (84) 
Ep 
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in terms of the on-sheh self-energy eigenvalue So(p). When compared with the non-relativistic HF 



energy (see eq. ( |207] ) in Section 4.2) it appears that, beyond the different expressions of the self- 
energy functions that hold in the relativistic and non-relativistic frameworks, an extra multiplicative 
factor m/E-p occurs in the relativistic case. 

Once the HF energy e(p) is known to first order in /^/m^, we expand as well the renormalized 
spinor s, namely 



/m(p) 
^(P) 



u(p,m(p)) 



E(p) + m(p) 



X 
a ■ p 



(85) 



^(p)+m(p) 

Actually, for later use in the calculation of the hadronic tensor, it is convenient to expand the 
spinor multiplied by the factor . / . 

V ^(p) 

Thus we start by expanding the Dirac mass in eq. (|75| ) around the on-shell energy, obtaining 
m(p) = m[l + Ao{p) + Co(p)] + 0{f/ml) ~ m + So(p) , (86) 
where use has been made of the on-shell self-energy in eq. (^2|). Likewise, using the HF equation 



(eq. (|6^)), the Dirac energy E{p) defined in eq. (|7^) is given by 
1 - B 



E{p) 



l-C 



e(p):^ [l-i?o(p) + Co(p)] 



m 



i?p + — So(p) 



m 



Ep + ;^So(p) ^ e(p) 



(87) 



After some algebra the following first-order expressions are obtained 



'E + m 
2E 
1 

E + fh 



I m + Ep I Ep-mT,o 

2Er. - 



1 



m + Er, 



1 



+ 
So' 

Er. 



2Ep Ep 



(89) 



Inserting eqs. (|88| ) and (^) into the renormalized spinor in eq. (|85| ) we get 



/m(p) 
^(P) 



u(p,m(p)) 



lm + E^ 



2Er. 



1 + 



Ep-m T^Q 

2En Er^ 



• P 

m + Ep 



X 

\ E„ , 



X 



I m Eo rm m + Ep 



Er. 



m 



2^. 



-X 



Ep + m cr ■ p 

2E„ m + E^ 



(90) 
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Since 



m)x 



{Ep-fo - m) 



X 

o- • P 

m + E-. 



-X 



(91) 



the first-order (in f"^ /m^) renormalized spinor can be cast in the form 



^u(p,m(p)):.^/- 



So(p)£;p7o-m 
^i(p) + — ^ ^i(p) 



2^. 



(92) 



The above expansion transparently displays the effect of the self-energy on the free spinor m(p). In- 



deed the second term in the square brackets of eq. (92) corresponds to a negative-energy component 
with momentum p. In fact, the Dirac equation for a positive-energy spinor is given by 

(p • 7 -h m)u(p) = Sp7ou(p), with Ep>Q . (93) 

Now if we apply the operator (p • 7 -|- m) to the spinor (-E'p7o — m)ti(p) we obtain 

(p • 7 -Fm)(£^p7o - m)ii(p) = p • 7(£;p7o - ■m)u(p) m(£'p7o - m)u(p) 

= (-^p7o - m)p • 7u.(p) + m(£'p7o - m)u(p) 

= (-£^p7o - m)(£;p7o - m)u(p) + m{Ep-fo - m)u{p) 

= -EpjoiEpjo-m)u{p). (94) 

Hence {Ep-fo — m)u is an eigenvector of the free Dirac Hamiltonian with eigenvalue —Ep. Therefore 
the operator -E'p7o — m transforms a positive-energy spinor u{p) into a negative-energy one. 

Moreover, it is useful to write down the correction to the free spinor (see eq. (^)) in an 
alternative form. Using the identity in eq. ( p^ ) we can write 

- m){Epjo - m)u{p) = 2Ep-foiEpjo - m)n(p) . (95) 

Multiplying by [2Ep{P — m)]"^ we then obtain 



Epjo - m 

2Er. 



n(p) 



1 



P — m 



7o(£^p7o - m)u{p) . 



Hence the second term in the square brackets of the right-hand side of eq. 
form 



(96) 

can be recast in the 



Sq £'p7o - rn 
Ep 2Ep 



u{p) 



1 



Ep P-m 



7o(£^p7o - m)u{p) 



1 



m 



f -m 



1 - —70 S(p)«(p) 



Sf{p) (1 - ^70 1 S(p)n(p) . 



(97) 
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The first term in eq. (|97D, /S'i?(p)$](p)n(p), corresponds to the one that enters in the original 
(divergent) self-energy current for a nucleon on-shell (eq. (p9D). The subtracted term, with the 



factor -^70 inserted between the propagator and the self-energy, cancels the divergence and yields 
a finite result. Thus it can be viewed as a "recipe" to renormalize the self-energy current. 

We turn now to an expansion of the field-strength renormalization function defined in eq. (66). 
For this purpose we use eqs. (|86| , |87|) , obtaining 

^2 ^ p2 Q(J 



^2(P) 



l + 5o(p) + 



m 



Ep dpo 



dpo Ep dpo 



J Po=Ep 



which implies that 



Z2(p) ~ 1 + -a(p) 



with 



a(p) = Bo(p) + 



dA 
Ep dpo 



dB p2 dC 



dpo Ep dpo 



(98) 



(99) 



(100) 



Po=Ep 



Hence, collecting the above results and inserting them into eq. (^), we get to first order 



/m(p)^ 
R)(p) 



m 



/ N So Ep-fQ -ml r \ r \ 
u{p) + — i^— n(p) + -a(p)u(p) 



Er. 



2Er. 



(101) 

Jp ^l^jp 

Thus, within the OPE approach the renormalized HF spinors in the nuclear medium are charac- 
terized by two new elements with respect to the bare m(p): the term ^ti(p) introduces 
negative-energy components in the wave function, while the term ^a(p)'u(p) comes from the field- 
strenght renormalization which modifies the occupation number of the single-particle states. 

Using the above expressions for the renormalized spinors, we now expand the renormalized 
one-body current matrix element to first order in /^/m^, getting 

' a(h). 



/m(p) /m(h) 



E{p) V E{h) 



i^F(P'h) ^ 



m m 



Ep Eh 



u(p) 



Eh 2£'h 



^ So(p)^p7o-"^ p^ ^ a(p)p^ 



Er, 



2Er, 



u{h) 



m 



(102) 



\/EpEh 

In eq. ( 102 ) the term jp^ represents the usual one-body current matrix element evaluated with 
free spinors, i.e., 



j^B(p,h) = n(p)r^(Q)n(h) , 



(103) 
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whereas j^g^ is a new renormalized self-energy current matrix element that includes the effects of 
the renormalization of the spinors. It can be decomposed according to 

ii?Si?(P. h) = Ji^5Ei(P' h) + JfiS£2(P> h) , (104) 
where j^^Ei embodies the correction arising from the new spinor solution of the Dirac equation in 
the medium and j'^gE2 stemming from the field-strength renormalization function \/^ in 

the medium. Their expressions are the following: 

f^SE^i^) = n(p) [r^ ^°^^^ " + ^°^P^ " n(h) (105) 

"a(h) +a(p)] 

JosvP'h) . (106) 



2.5 Gauge invariance of the theory 

A crucial feature of the present theory is that the hadronic tensor, computed either through the 
p-h matrix elements or through the polarization propagator, is gauge invariant. This may be 
somewhat surprising because, as shown in Appendix B (see also |jl9|), current conservation is 
already obtained at the level of the MEC and correlation p-h matrix elements: hence the one-body 
current p-h matrix element also has to be independently conserved. This however occurs only in 
zeroth order of perturbation theory. To be dealt with properly, the situation clearly requires the 
renormalization of the p-h energies and of the Dirac spinors (see previous sections) . Only then does 
it become possible to set up a renormalized SE current which leads to a hadronic tensor coinciding 
with the one obtained through the polarization propag ator 111. 



As shown in the previous section, the renormalized HF current matrix element, expanded to 
first order in /^/m^, has been split into the usual one-body current and into a new renormalized 
self-energy current. In order to be consistent with the one-pion-exchange model, we should add 
the contribution of the seagull, pion-in-flight and vertex correlation currents corresponding to the 
diagrams shown in Fig. 2(a-e). We point out once more that the self-energy diagrams (f) and (g), of 
Fig. 2, corresponding to insertions in external legs, should be disregarded in computing amplitudes 
(or currents) in perturbation theory. Rather, their contributions should be taken into account via 
renormalized energies and spinors as solutions of the relativistic HF equations. We have expressed 
the latter, to first order in /^/m^, in the form of a new current operator (denoted as RSE current). 

Then the total current in our model reads 

j^(p,h) = jg^(p,h) +i^p^(p,h) , (107) 
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where Jope embodies the seagull, pion- in- flight, vertex correlation and renormalized self-energy 
currents, namely 

j&PE = + j!; + ji^c + Jrse ■ (108) 
In what follows we shall prove the gauge invariance of this current to first order in /^/m^. In 



so-doing it is crucial to take into account not only the full current in eqs. (107, lOq ), but also the 
first-order correction to the energy of the particles and holes due to the self-energy interaction in 
eq. (|8^ . In other words, for a given momentum transfer q = p — h, the energy transfer should be 
computed as the difference between the particle and hole HF energies and not using the free values 
Ep and E^. Thus the energy transfer is 

TTl TTl 

UJHF = Ep-E^ + — So(p) - ^So(h) (109) 

and the associated four-momentum transfer is Q^p = {ujhf,^)- To make the following discussion 
clearer we denote with Qhf the HF four-momentum and with cohf the HF energy transfer, to 
distinguish them from the on-shell values Q and uj. 

Divergence of the one-body current 

The divergence of the zeroth-order one-body current computed using the HF four-momentum 
transfer Qhf is given by 

QHF,^j^BiP, h) = u{p)QHF,^r'' (Qhf) n(h) = u{p)Fi (Qhf) ^HFu{h) , (110) 

where the nucleon vertex T^(Qhf) is also evaluated at the momentum transfer Qhf- Because of 
u(p)Qu{h.) = 0, only the first-order contribution arising from the self-energy correction survives, 
namely 



Qi?F,/.JoB (P,h) = u{p)Fi{Q) 



I^So(p) - I^So(h) 



7on(h) . (Ill) 



In the above the Dirac form factor Fi is computed at the unperturbed value Q'^, since we disregard 
second-order contributions. Note that the one-body current itself is not gauge invariant — its 
divergence yields a first-order term which turns out to be essential for the gauge invariance of the 
full current, as we shall see below. 
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Divergence of the MEC 

The seagull and pionic Ip—lh currents given in eqs. ( [25|p6| ) are already of first order in f'^/m'^; 
thus in computing their divergence we use the unperturbed value of the energy transfer, neglecting 
a term of order O Using the free Dirac equation and exploiting the kinematics we obtain 



Q,jT(p,h) 

f -/X ^ m (2(K ■ P -m ^) 2(K ■ H - m ^)] , 

= -^^1 -3..^(P)r.r. E ^ [(^p.Ky.ml - {K-HY-ml\ ^^^^ 



(112) 



P r.v- , ^ V- f 2m(^-m) 2m{^ - m) \ , 



k<fc 



(113) 



(114) 



In deriving these equations we have used the relations Q^{Q + 2H — 2K)^ = —2K ■ Q and 

1 1 _ -2P-Q 

{K - Hf - ml ~ {P - Kf - ml ~ [{K - Hf - ml] [(P - - ml] ' 

Upon addition of Eqs. ( |112| , |113| ) the terms containing ^ cancel, leaving for the total divergence of 
the seagull and pion-in-flight the expression 

f V s V- m{2iK-H-m'^) 2(K-P-m'^)\ , , ^ 

which can be further simplified by exploiting the self-energy of eq. ( |4^ ) for on-shell momenta. One 
finally obtains 

Q^^{j'^ + = '-FYe3abuip)ran[np) - S(h)]tx(h) . (116) 

Divergence of the vertex correlation current 

Starting from the Ip-lh matrix element of the VC current in eq. (^) and applying the Dirac 
equation, we get 

f 2 I K + m 



F 



f 2 1 X -\- m 

^^(p)r.P... E^^75(;r- /^) (^_%_^. 75(/^- PHI.) , 

(117) 
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where we recognize again the expression of the self-energy matrix in eq. (|4l| ) . Since the Dirac form 
factor can be spHt into an isoscalar and an isovector component according to 
1 



(118) 



which yields 

TaFiTa = 3Fi + e^abTaU , (119) 

the divergence of the VC current written in terms of the self-energy function reads 

(120) 



Q/.j{^c(P' h) = n(p) ( Fi + -F^e^abTan ] [S(h) - S(p)] u{\i) . 



Comparing this result with eq. ( |116D we note that the term above containing t^ah'^a'^h cancels with 
the MEC contribution. Hence 

bMEc(P' h) + j(;c(P' h)] = n(p)Fi [S(h) - S(p)] u{h) . (121) 

The above relation just expresses the Ward-Takahashi identity |46| relating the full vertex correc- 
tion, namely MEC plus VC (diagrams 2 (a)-(e)), to the self-energy matrix element. 

Divergence of the RSE current 

Finally we compute the divergence of the renormalized self-energy (RSE) current defined in 
eqs. (104), (|105D and ( |106D . For this purpose we first note that the divergence of jjisE2 vanishes to 
first order because it is proportional to the OB current. Hence we write 



T.o{h) Ehjo - m So(p)£'p7o-m 



Eh 2Eh 
Using the relation u(p) ^u{h) = and 



+ 



2Er. 



u(h). 



(122) 



u{p) @7ou(h) = «(p)2(m7o - EMh) (123) 

w(p)7o ^u{h) = u{p)2{Ep - m7o)n(h) (124) 

it is straightforward to obtain 

Q^j^S^(p,h) = n(p)Fi [S(p) - S(h)] u{h) +u{p)Fi 

Remarkably the first term of this equation cancels with the divergence of the MEC plus the VC 
current, given by the Ward-Takahashi identity in eq. (|121|), whereas the second term cancels with 



m 



m 



— So(h) - — So(p 



u(h) . (125) 



the divergence of the OB current in eq. (111). We have thus proven that, within the present model 
up to first order in /^/m^, the total current in eq. (107) satisfies the continuity equation, namely 

QheAJob + Jmec + JVC + Jhse) = • (126) 
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(a) (6) (c) (rf) 




(e) (/) (9) [h) 



Figure 6: Feynman diagrams of the free (a) and first-order pion-in-fiight (b), seagull (c and d), 
vertex correlation (e and f) and self-energy (g and li) polarization propagator. 

2.6 Nuclear hadronic tensor and electromagnetic response functions 

In this section we compute the electromagnetic inclusive response functions for one-particle emission 
reactions within the RFG model. As discussed in previous sections, the p-h matrix elements 
corresponding to the different pionic diagrams are all well-defined except for the self-energy term 
which diverges, and consequently needs to be renormalized. In what follows, we evaluate the 
hadronic tensor starting from the current p-h matrix elements in the case of the one-body, MEC 
and vertex correlation diagrams. These are shown diagrammatically in Fig. 6. On the contrary, 
for the self-energy diagrams we calculate the hadronic tensor in two at first sight different ways: 
on the one hand, from the polarization propagator 11^^ (see Appendix C), and on the other, using 
the renormalized, well defined, SE p-h matrix elements (Appendix D). We prove that the two 
formalisms are equivalent. 
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The formalism of the nuclear hadronic tensor set up with the p-h matrix elements has been 



presented in detail, within the RFG model, in previous papers [24, Hence, here we simply sum- 
marize the results needed for later discussions. Before starting the analysis of pionic contributions, 
we recall the analytic expressions for the OB, leading order electromagnetic responses of the RFG 
(see, for example, |^, ^] for details): 



(127) 



where p and n refer to protons and neutrons, respectively, and, for Z = N, 
3Z 



Ro{q,uj) 



AmKrjp 



{ef - eo)9{eF - eo) 



with 



eo = Max je^ - 2A,Ky 1 + - - A| . 
In the above the usual dimensionless variables 



A = , T = ^. 

2m 4m^ 



n = VF = — , 

Zm m 



Ep 

m 



have been introduced and Ep = \jkp + m? is the Fermi energy. The functions U ' in eq. 
are 



(128) 
(129) 
(130) 



A 



A 



1 + r 



where 



A 



e% + epEQ + ^o) + {^F + eo) + A' 



(1+r) 



(131) 
(132) 

(133) 



2.6.1 MEC and vertex pionic contributions 

The hadronic tensor that arises from the interference of the single- nucleon, OB current, Jq^, with 
the one-pion-exchange current j^, with a = s (seagull), p (pion- in- flight) and vc (vertex correlation), 
is for the RFG model with Z = N (see eq. (|193|) below) 



3Z 



pq Jho 



kp /•2tt 

hdh{uj + ^h) / #h 



m 







EpEh 



2Re[i^^(p,h)*j:(p,h)] 



(134) 
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where jo^(p, h) = M(p)r'^M(h) is the smgle-nucleon p-h matrix element with the electromagnetic 
nucleon current from eq. ( p2[ ) and j^(p, h) is the p-h matrix element for the seagull, pion-in-flight 
or vertex current as given in eqs. (|25|), (p^) and (pg|-]30|), respectively. 

Note that in eq. ( |134| ) the integral over the hole polar angle, cos 9h, has been performed explicitly 
by exploiting the energy-conserving 5-function. This fixes the minimum momentum of the hole 
according to 



ho = my Eq — 1 . (135) 
Moreover, the hole three-momentum 

h = h (sin 60 cos (j)^, sin 60 sin cp^, cos 60) , (136) 

involved in the hadronic tensor, must be evaluated for the following specific value of the polar angle 

cos0o = ^^^, (137) 
r]K 

with 7] = h/m. 

The hadronic tensor, as was the case for the current, can be also split into isoscalar and isovector 
parts, since there is no interference between the two isospin channels. 

An important issue relates to the form factor of the ttNN vertex, F^, which incorporates some 
aspects of the short-range physics affecting the pionic correlations. In all of the above expressions 
has not been explicitly indicated for sake of simplicity. In the analysis of the gauge invariance 
at the level of the particle-hole channel, performed by deriving the contribution to the continuity 
equation of the isoscalar and isovector SE, VC and MEC p-h matrix elements, is presented. There, 
it is shown that the SE and VC contributions cancel in the isoscalar channel, in contrast with the 



non-relativistic result [24|, where the SE is by itself gauge invariant. Furthermore, the SE and VC 
contribution in the isovector channel is exactly canceled by that of the MEC (seagull and pion-in- 
flight). It is crucial to recall that the inclusion of F^r in the p-h current matrix elements is not without 
consequences in connection with gauge invariance. In fact, in this case, the model is not gauge 



invariant unless new terms are added to the MEC (see ||5l| , |52| , |53| , p^ , |55| and |5(:, ^ for recent 
work on the restoration of current conservation in model calculations). Lacking a fundamental 
theory for F^r, in the calculations reported in this work we use the phenomenological expression 



2 
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with A = 1.3 GeV. As long as the dependence upon A is not too strong the gauge invariance of the 
theory should not be too badly affected. Within a non-relativistic approach for the pion currents, a 
detailed discussion on the breakdown of the gauge invariance induced by Ftt, and on the dependence 



of the responses upon the cutoff value can be found in |55|| . 



In [58, 59] the effects of the MEG upon the transverse response in a non relativistic shell model 
for finite nuclei were studied as a function of the cutoff A. 

2.6.2 Relativistic self-energy responses 

As already discussed in previous sections, a crucial point to be emphasized is that the self-energy 
p-h matrix element, eq. (pg|), is divergent. Hence it cannot be used directly in the evaluation of the 
hadronic tensor. Instead one should use renormalized spinors with the corresponding renormalized 
energies. Above we have taken account of the effect of renormalization to first order in /^/m^ by 



introducing an extra term in the current: the RSE current defined in eq. (104). In addition there 
is also a 0(/^/m^) modification of the energy of the particles, eq. (p4[). These two modifications of 
the free current and energy in turn give a contribution to the hadronic tensor of order 0(/^/m^), 
which we will refer to as renormalized- self energy contribution (RSE), which is of the same order 
as the MEC and VC currents and should be included in any consistent calculation to first order in 
/■^/m^. In addition this contribution is needed for the gauge invariance of the results. 

In what follows we derive the RSE contribution to the nuclear response functions. This RSE 
contribution should replace the SE Feynman diagrams shown in fig. 6 (g),(h). As a matter of fact, 
these two diagrams can be computed using the polarization propagator formalism (see Appendix 
C and Ref. [^]), where one does not need to appeal to renormalization since the SE diagrams 
are finite in this case. Our goal is to show that the results for the response functions obtained in 
the two ways coincide, although they stem from different approaches. This is proved in Appendix 
D. The RSE contribution, therefore, can be identified with the contribution coming from the two 
diagrams (g), (h) of fig. 6. 

The one-body hadronic tensor in HE approximation reads 

X (^(cj + e(h) - e(p))^(fcF - h) , (139) 

where p = h + q and j^p.(p,h) is the one-body HE current in eq. (^) computed using the 
renormalized HE spinors and HE energies of the particle and the hole. 
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Next we use the expansions in eqs. (102) for the current jup and ( p4D for the HF energies. In 
addition we expand the energy delta function to first order in /^/m^ according to 

d5{u; + Eh- Ep) 



6{uj + e(h) - e(p)) ~ 6{uj + Ei,- Ep) + 



duj 



m 



m 



So(h) - ;^So(p) 



(140) 



Inserting all of these relations into the hadronic tensor and neglecting terms of second order we get 
for the diagonal elements of the hadronic tensor ^ 



Wi^^pico, q) ^ W^^siiO, q) + AH^j^^^(u;, q) 



(141) 



(the summation convention is not in force in eq. (|141| )), where Wq^{uj, q) is the usual OB hadronic 
tensor of a RFG, i.e., 



and AWj^g^{uj,c[) is the first-order self-energy correction 



(142) 



m 



+ \f0BiPM\' 



So(h) - ^Eo(p 



2Re j^^ip, h)*j^5g(p, h)5{uj + E^- Ep) 
d 



Eh 



Er, 



dio 



6{oj + - Ep) } e{kF - h) . (143) 



In eq. (|143| ) the first term corresponds to the interference between the OB and the RSE currents, 
while the second one, which shifts the allowed kinematical region because of the derivative of the 
energy delta function, is due to the modification of the nucleon energies in the medium. 
Carrying out the spin traces for the single-nucleon current 

1 



SpSh 



Am? 



Tr {T^'{Q){P + m)T^'{-Q){P + m)} 



(144) 



we get for the renormalized self-energy response function 



= V 



d^h 



(27r)3 AEpEh 



Tr r^(Q) 



Eo(h) Eh7o - ^ a(h) 
E\i 2£'h 2 



{P + m)T^'{-Q){P + m) 



+Tf'{Q){P + m)r'^(-Q)(/' + m 
+r''{Q){p + m) 



T.o{h) Sh7o -m ^ a(h 



Eh 2£'h 



So(p) Epjo - m ^ a{p) 
Ep 2Ep 2 

rn-Q){P + m) 



^We only work out the diagonal elements of the hadronic tensor, since these are the ones that contribute to the 
unpolarized inclusive longitudinal and transverse response functions. 
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So(p) -Ep7o - m ^ a(p) 



x6{uj + Eh - Ep)e{kF - h) 
+ y j j^,j^T,{T^{Q){P + m)T^{-Q){P + m)} 



'J + m) 



d^h 1_ 



^^So(h) - — So(p) 1 —6{uj + Eh- Ep)e{kp - h) . (145) 



More precisely, one should add two copies of eq. ( 145 ), one with the form factors appropriate to 
the proton and one to the neutron. 

In Appendix D we show that this contribution to the response function is identical to the 
one obtained in Appendix C by computing the imaginary part of the polarization propagator 
corresponding to the two SE diagrams (g), (h) of Fig. 6. This identity is not trivial: indeed in the 
case of the polarization propagator the response functions, with the Fock self-energy dressing the 
particle and the hole lines, are computed by representing the product of two nucleon propagators 
as the derivative of a single one to deal with the presence of a double pole in the integrand. In 
the present paper the problem has been solved differently. First the entire perturbative series with 
Fock self-energy insertions has been summed up and then the result has been expanded to first 
order, thus obtaining a finite first-order current operator. Because of the equivalence of these two 
procedures we are confident about the validity of the results we have obtained for the self-energy 
contribution to the nuclear responses. 

2.7 Analysis of results 

In this section we report the numerical results obtained for the pionic MEC (pion-in-fiight and 
seagull) and for the correlation (vertex and self-energy) contributions to the quasielastic peak 
(QEP) in the Ip-lh sector. The calculation is fully-relativistic. We have taken Z = N = and 
set kp = 237 MeV/c, which is representative of nuclei in the vicinity of ^'^Ca. 



The 5-dimensional integrations of the MEC and correlation responses implicit in eq. (134) have 
been performed numerically. The reliability of the numerical procedure has been proven by checking 
that the free RFG responses coincide with their analytic expressions (see, e.^., []35| ). 

2.7.1 MEC 

We start by analyzing the effects introduced by the MEC. These are presented in Figs. 7 and 8 
where we show the longitudinal (left panels) and transverse (right panels) response functions versus 
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Figure 7: Longitudinal (left panels) and transverse (right panels) electromagnetic response func- 
tions versus u. Dashed: free RFG; solid: RFG+MEC contribution. Here and in all the figures 
that follow, unless explicitly indicated, the nucleus is ^°Ca, corresponding to a Fermi momentum 
kp = 237 MeV/c. 
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Figure 8: Separate MEG contribution to the longitudinal (left panels) and transverse (right panels) 
responses. Dashed: pion-in-flight; short-dashed: seagull and solid: MEG (pion-in-flight + seagull) 
contribution. 
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the transferred energy lo for four different values of the transferred momentum q: 0.5, 1, 2 and 3 
GeV/c. First, in Fig. 7 we compare the free RFG responses (dashed) with the responses obtained 
including the global MEG contribution (solid). As shown, while for the longitudinal responses 
the MEG are hardly visible, in the T channel they contribute somewhat more, typically by about 
5-10%, depending upon q and oj (see discussion later). 

In Fig. 8 we display the separate pion-in-flight (dashed) and seagull (short-dashed) contribution 
to and ii^ for various values of q. The total MEG (seagull + pion-in-flight) contribution is 
also shown (solid line). In the transverse channel (right panels) it appears that the seagull term 
is always larger than the pion-in-flight term, a dominance that increases with q and reflecting the 
spin nature of the photon-MEG interaction. Moreover, whereas the pion-in-flight term is always 
negative, the seagull changes sign with w, inducing a (mild) softening of the response, i.e., a shift to 
lower energy. Within the longitudinal channel (left panels), the seagull term, now always negative, 
also dominates. Note however that the relative difference between the seagull and pion-in-flight 
contributions is not as large as in the previous case. Moreover, the behavior of the seagull and 
pion-in-flight terms in the longitudinal channel as q increases displays a different pattern from the 
one shown in the transverse channel, since for high q the pionic current is not negligible compared 
with the seagull one. 

To complete this discussion we briefly comment on the MEG dependence upon the momentum 
transfer q and the Fermi momentum kp, associated with scaling of first and second kind, respectively 



(see [|q, 61, Bl). 



In [19 1 we have explored in detail the evolution with q of the MEG in the transverse channel (as 
they are negligible in the longitudinal channel). We have proven that their relative contribution 
to decreases with q, but does not vanish for large values of q. In fact, the relative MEG 
contribution decreases in going from 0.5 to 1 GeV/c, but then it rapidly saturates at or slightly 
above q=l GeV/c, where its value stabilizes, typically around 10%. Thus, one can conclude that at 
momentum transfers above 1 GeV/c scaling of the first kind is satisfied for the MEG contributions 
considered in this work. Moreover, for high q the MEG almost vanish for w in the vicinity of the 
QEP. 

A detailed analysis of the kp dependence of the MEG contribution in the transverse response 
has also been presented in |19]. The MEG contribution is found to grow with fc^?, in contrast with 
the free response which decreases as k'^^ . It is also shown that the two-body MEG processes violate 



40 



the second- kind scaling by roughly three powers of kp. This effect is a rapid function of the Fermi 
momentum (or equivalently, of the density): for example, if one considers the cases ^H/^He/heavy 
nuclei with Fermi momenta of approximately 55/200/260 MeV/c, respectively, then the Ip-lh MFC 
contributions amount to 0.1/5/10% of the total transverse response, respectively (normalizing to 
10% for the heavy nucleus case). 

2.7.2 Correlations 

In Fig. 9 we display the vertex correlation contribution to the longitudinal and transverse responses 
by comparing the free RFG responses (dashed) with the responses obtained including the VC 
contribution (solid). As noted, the VC action, while substantial in both the longitudinal and 
transverse channel, is actually dominant in the former by roughly a factor of three. This outcome 
relates to the minor role played by the isoscalar contribution in the transverse response, in turn 
due to the smallness of the isoscalar magnetic moment. 

The evolution with q of the VC in the longitudinal and transverse channels has been discussed at 



length in |1£]. Let us summarize the basic findings. First, the VC do not saturate quite as rapidly 
as the MFC, although their behavior is rather similar and saturation again occurs somewhere above 
q = 1-1.5 GeV/c: thus, once more, scaling of the first kind is achieved at high momentum transfers 
for these contributions. Moreover, similarly to the MFC case, for high q the VC almost vanish 
around the QFP. 

Finally, the vertex correlations are found to grow with kp, much as the MFC do. From a 
semi-relativistic point of view, we find a behavior that goes as kp. The basic conclusion is similar 
to that made above for the seagull contribution and hence for the total MFC at high q, namely, 
scaling of the second kind is badly broken by effects that go roughly as kp. 

The role played by the SF contribution is displayed in Figs. 10 and 11. In Fig. 10 we show the 
total RFG + SF responses (solid line) compared with the free RFG responses (dashed). Note that, 
in contrast with the MFC and vertex correlations, which mostly contribute to only one channel 
(transverse and longitudinal, respectively), the impact of the self-energy on and is similar, 
leading in both cases to a softening of the responses for high q. 

The separate particle (dashed) and hole (short-dashed) SF contributions to the longitudinal 
and transverse responses are presented in Fig. 11. Here, also the total SF contribution (solid) is 
displayed. We observe that the self-energy contribution results from a quite delicate cancellation 
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Figure 9: Same as Fig. 7 but for the vertex correlation. Dashed: RFG responses; sohd: RFG + 
VC contribution. 
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Figure 10: Same as Fig. 7 but for the self-energy. Dashed: RFG responses; solid: RFG + SE 
contribution. 
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Figure 11: Particle (dashed) and hole (short-dashed) contributions to the longitudinal (left panels) 
and transverse (right panels) self-energy. The solid line represents the total SE contribution. 
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between the responses having only the particle or only the hole dressed (Fig. 11). This was already 
pointed out in |51] within the framework of a treatment in which relativistic effects were partially 
incorporated and it is now confirmed within a fully-relativistic context. 

Whereas this cancellation is very substantial at q=0.5 GeV/c, as the momentum transfer in- 
creases the imbalance between the two contributions grows. Indeed the response associated with 
the particle self-energy is suppressed by the form factors and by the pion propagator, but that 
coming from the hole self-energy is not. As a result, for q > 2 GeV/c the total self-energy response 
is almost entirely due to the hole dressing and induces a moderate softening to the free response. 
Note that the SE contribution does not vanish on the borders of the response region. Moreover 
for high values of uj (close to the upper border) it becomes very large (Fig. 11) and yields a sig- 
nificant lowering of the upper tu limit in the responses. This clearly points to the insufficiency 
of a ffist-order perturbative treatment in this kinematical region, an effect already present in the 



partially relativized analysis of [51| and emphasized by our fully-relativistic calculation. Therefore 
the summation of the full Fock series becomes necessary near the upper boundary of the response. 
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Figure 12: The on-shell self-energy Eo(p)/-Ep defined in eq. ( p2| ) (solid line) and the field-strength 
renormalization function a{p) given in eq. (|100|) (dashed line) plotted versus the momentum p. 



The analysis of the scaling and superscaling properties of the self-energy correlations has been 
presented in [|l9| ]. In accord with the above, the particle contribution decreases with q, going to 
zero at g ~ 2 GeV/c, whereas the hole contribution, although also decreasing with q when not 
too high, saturates for q >1 GeV/c (see Fig. 11). As a result the total self-energy grows with q 
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and transverse (right panels) responses plotted versus uj. The separate contributions of the current 
^RSEi ^^"^ particle (solid) and hole (dotted) and of the current jjisE2 particle (dashed) 

and hole (dot-dashed) are displayed. 
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Figure 14: The contributions of the first (solid) and second (dashed) term in eq. ( |143D to 
longitudinal (left panels) and transverse (right panels) responses. 
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in the range g=0.5-2 GeV/c, then stabihzes typically at about 30-40% of the free response to the 
left of the QEP, thus inducing an important softening of the longitudinal and transverse responses. 
In summary, again scaling of the first kind is achieved at momentum transfers somewhat below 2 
GeV/c. Finally, we also prove that the self-energy relative contribution grows with kp, although 
not uniformly in the scaling variable i/j (see pO] , |6l| , |6^ ) — recall that in the first-order analysis 
presented in this paper the edges of the response region are not treated adequately for the self- 
energy contribution and thus should not be taken too seriously. Where the self-energy contribution 
is correctly modeled (away from the edges) we again see breaking of second-kind scaling by roughly 
kp. 

In what follows we explore the impact on the responses of the new currents jfisEi jRSE2 
that arise from the enhancement of the lower components of the spinors and from the field strength 
renormalization ■^/^(p), respectively. In Fig. 12 we show the on-shell self-energy (solid curve) and 
the field strength renormalization function (dashed curve) given by eqs. (|8^ ) and ( |100| ), respectively. 
The explicit expressions for 5]o(p) and a{p) are derived in Appendix E. The So(p) obtained here is 
in good agreement with the results of and its effect on the single-particle energy in eq. ( [8^ ) and 
on the effective mass in eq. ( ^6|) is very small (less than ~ 3%). Note that q, which is linked to the 
current j'^gE2 °^ ( |f06| )i is much smaller than So(p)/i?p, which enters in jj^g^i through eq. ( |105D . 
Thus the effect of the enhancement of the lower components of the spinors dominates over the field- 
strength renormalization. This is very clearly seen in Fig. 13, where the various contributions to 
the longitudinal and transverse responses stemming from jj^gpi and jjigp2 displayed versus the 
transferred energy u for momentum transfer q = 0.5, 1, 2 and 3 GeV/c. It is evident that the effect 
of jfisE2 is negligible with respect to that of jfisEi- The separate contributions of the particle 
and hole self-energies are also shown: as q increases the contribution of the particle is suppressed, 
whereas the one of the hole survives. 

In Fig. 14 we compare the contribution to the longitudinal and transverse responses due to 
renormalization of the wave functions (solid) with that arising from renormalization of the energies 
(dashed). The effect linked to modification of the energy due to the medium is the dominant one, 
the other being very small, especially for large values of q. 

To complete this section we display in Fig. 15 the separate contributions of seagull (dashed), 
pion-in-flight (dot-dashed), VC (dotted) and SE (solid) to the longitudinal and transverse responses. 
Worth pointing out is the oscillatory behavior versus uj of the vertex correlations, which induces 
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Figure 15: Separate pion-in-flight (dot-dashed), seagull (dashed), vertex correlation (dotted) and 
self-energy (solid) contributions to the longitudinal (left panels) and transverse (right panels) re- 
sponses. 
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Figure 16: Longitudinal and transverse vertex correlation responses versus uj in the pseudovector 
(solid) and pseudoscalar (dashed) -k-N coupling. 
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Figure 17: Longitudinal and transverse responses versus lo including all first-order contributions 
(solid) compared with the free result (dashed). 
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a hardening of the responses. In addition the seagull and vertex correlations tend to cancel in the 
transverse channel, especially for low values of q, whereas for higher q the MEC dominate. Note 
that in the T channel both the seagull and VC exactly vanish at the same value of to, the latter 
coinciding with the QEP for high momentum transfers, as said above. It is also important to point 
out that the net effect introduced by the SE contribution is in general the largest one for transfer 
momentum values q > 1 GeV/c. Within the L channel, the pionic correlations (VC and SE) clearly 
dominate over the MEC. In the transverse channel, apart from the SE contribution which seems 
to dominate for high q, the seagull term is clearly more important than the VC one, whereas the 
pion-in-flight only enters for q not very high. 

Up to now we have considered a pseudovector coupling for the pion, eq. ([T^. We now shortly 
investigate the effects on the responses of using a pseudoscalar pion-nucleon Hamiltonian 

'^i^NN = W'^l^^'aTalp (146) 

instead of the pseudovector one. For on-shell nucleons the Hamiltonians in eqs. (|l^) and ( |146D are 
equivalent provided f/m-,^ = g/{2m), but for off-shell nucleons this is not so. Among the diagrams 
considered in our approach the only one involving off-shell nucleons is the one associated with the 
vertex correlations (Figs. 6e and f). Hence in Fig. 16 we compare the VC contribution to and 

obtained with the pseudovector (solid) and pseudoscalar (dashed) couplings. The difference 
between the two is especially sizable in the transverse channel (where the impact of VC is smaller) 
and increases with the momentum transfer. 

In conclusion, in Fig. 17 we display the total responses in first order of perturbation theory and 
compare them with the zeroth-order ones (free responses) for several momentum transfers. Here one 
assesses the impact of the global two-body current contribution to the responses. First the overall 
effect of the two-body currents appears sufficiently modest to justify our first-order treatment. 
Next the softening at large q appears to be common to both L and T channels, whereas at low q 
the longitudinal response displays a hardening that is absent in the transverse one. Also evident is 
the already-noted nearly vanishing of the two-body correlation contribution at the peak of the free 
responses. Finally the unrealistic dominance of the self-energy contribution on the upper border 
is apparent. Summarizing, the impact of the different first-order contributions — MEC, vertex 
correlations and self-energy — to the total responses are all comparable in size in the transverse 
channel (in the longitudinal one the MEC are negligible), their relative contribution ranging from 
~5 to ~15% depending upon the kinematics and the Fermi momentum. 
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3 Parity- violating electron scattering 

In this section we deal with the parity-violating (PV) effects arising from the weak interaction 
between the electron and the nucleus. Such effects, which are negligible in unpolarized electron 
processes, can be brought to evidence by measuring the asymmetry associated with longitudinally 
polarized electrons having opposite helicities, namely 

A — ~ (147) 

In this case the purely electromagnetic cross sections cancel out and one is left with the interference 
between the electromagnetic and neutral weak currents, corresponding to the exchange of a photon 
and a Z^, respectively. 



An important motivation of parity- violating experiments (see, for example, [64| for a general 
review and ||3^ for the foundations of PV quasielastic scattering) is the measurement of the single- 
nucleon form factors, in particular the strange and axial ones: for this reason most experiments are 
presently being carried out on light nuclei, where the uncertainties associated with the nuclear model 
are minimized. Other motivations exist for such studies: specifically, as discussed in the following, 
the PV response functions display a different sensitivity to nuclear correlations compared with the 
parity-conserving ones: hence they could not only shed light on the part of the problem concerned 
with nucleon (and meson) structure, but also are being used as a test of nuclear models. In the 
present work we provide no details for the underlying formalism used in PV electron scattering — 



those discussions can be found in [64|. Our focus here is rather to place in context the expectations 



for PV electron scattering of what role the modeling discussed above plays. 
3.1 General formalism 

The cross section for scattering of a polarized electron with helicity h reads 



In eq. ( [I48D G is the Fermi constant, 77^,^, W^'^ are the leptonic and hadronic electromagnetic 



tensors defined in Section 2.1 and fj^j^y, W^'^ are the tensors arising from the 7 — interference. 



Here terms containing two weak currents have been neglected. The interference tensors read 

flt^u = {av - haA) (k^KI + K'^K, - g^,K ■ K') + (a^ - hav)ie^,p^KPK"' (149) 



53 




Figure 18: Feynman diagrams of the free (a) and first-order pion- in-flight (b), seagull (c and d), 
vertex correlation (e and f) and self-energy (g and h) PV polarization propagator. The external 
wavy and dashed lines represent a photon and a Z'^ boson, respectively. 

with aA = —1 and ay = 4sin^ 6w — l, Ow being the weak mixing angle, and 

w'''' = Y.T.(f\Jernmir{f\JUQ)\i)m+^-Ef) , (150) 

i f 

JemiQ) -^wniQ) being the nuclear electromagnetic and weak neutral currents, respectively. 

When the difference of cross sections corresponding to opposite electron helicities is taken, the 
electromagnetic term in eq. ( |148| ) cancels out and the resulting PV expression reads 
da \ if d(T^+^ do-(-) \ 



dVl'i^duj J py 2 \dQ,'^duj dft'^duj J 

= Aqo-m VLRj,y{q,uj) + vtRav{i^^) + VTRvA{<i^^) ; (151) 

where 

A = , (152) 
2v 27ra 
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aM is the Mott cross section in eq. @), the leptonic kinematical factors vl and vt are given by 
eqs. and 



VT' = tan^ J- ( — ) +tan2 ^ . (153) 



In terms of nuclear response functions the asymmetry in eq. ( |147D reads 



3.2 PV response functions 

The PV response functions appearing in eq. ( |151| ) are hnked to the interference hadronic tensor in 
eq. ( |150| ) by the following relations: 



q 



(155) 



R^viQ,^) = aA{w^^ + W^^) (156) 
Rva{q,^) = -iav{w^^-W^^) . (157) 

The subscript AV in the PV responses denotes interferences of axial-vector leptonic currents with 
vector hadronic currents, and the reverse for the subscript VA. 

Within the context of the RFG model the interference hadronic tensor is 

= ^T^ / ^^^(^ + / #h E ^^2Re [ji:^{p, h)*j-„(p, h)] , (158) 

where the electromagnetic current includes both the single nucleon one-body and the two-body 
(MEG and correlation) currents discussed in the previous section, i.e. = Jq^ ~^ Jmec ~^ 3cor- 
In this work we include in the weak neutral current only the one-body contribution (see Fig. 18), 
namely 



J-; Jp, h) = n(p) ( F^Y + ^^''"'Qp + GaI^i" 1 u{h) , (159) 



where the Pauli and Dirac form factors are 

A = Ge + tGm (IgQ) 

l + T 

= (161) 



55 



Thus we neglect the direct couphng of a to the pion (important clues for the understanding 
of the weak-neutral sector of the MEC should be found in the study of pion electroproduction on 
the nucleon, where a is exchanged with the nucleon. This topic has recently been investigated 
in H). 

Within the standard model at tree level the weak neutral form factors are linked to the electro- 
magnetic ones by the following relations (possible contributions from the strange quark are neglected 



— these can be included in a straightforward way 1 35 , Q and do not provide the primary focus of 
the present discussions): 



G 



Ep{n) 



GMp{n) — PyGMpin) + Pv^Mnip) 



G 



Ap{n) 



0AGAp{n) + P\GAn{p) 



where 



V 



1-4 sin^ 



(162) 
(163) 
(164) 

(165) 



The one-body contribution to the three PV responses can be evaluated analytically in RFG, yielding 
(see, for example p5[) 



L,Ti 



00 



RvAiq^oj) 



ayRoiq, u;) Up {q, uj) + {q, u) 



where Rq has been defined in eq. ( |128| ) and 

if ^ A 

r 

VT 



G 



Ep{n)GEp{n) + 

A 



G Ep{n)G Ep{n) + 'T G M p{n)G M p{n) 



^p{n)il^^) — '^'tG Mp{n)G Mp{n) + _|_ ^ G Ep{n)G Ep(n) + tG Mp{n)G Mp{n) 



2Wr(l + r)GAfp(„)G^p(„)(l + A) 



with A given by eq. (133) and 



^ (ep + eo) + A 



1 



(166) 
(167) 

(168) 

(169) 
(170) 

(171) 



K V 1 + r 

The two-body contributions involve instead multidimensional integrals, to be numerically evaluated. 



56 



3.3 Results 

In this section we analyze the PV response functions labeled L, T and T' and the associated 
asymmetry for various values of the momentum transfer. In ||6^ results for the PV responses in a 
relativized continuum shell model where presented in the impulse approximation for finite, closed 
shell nuclei. Two-body currents where not included in that calculation. In |51 , ^| a semi-relativistic 



analysis of the PV responses has been presented, showing the dominance of pionic correlations in 
the longitudinal channel. Here we perform a fully-relativistic calculation, which confirms the above 
findings, and extend them to higher values of the momentum transfer. 

In Fig. 19 we display the PV responses for four values of the momentum transfer q. The dashed 
line corresponds to the free RFG, the dotted line includes MEG and vertex correlations, whereas 
the solid line also includes the self-energy contribution. 

One observes that in the T and T' channels (central and right columns) the main effect arises 
from the self-energy, which tends to soften the response function, similarly to what happens in 
the parity-conserving case, whereas the MEC and VC effect is very tiny. Note also that the 
axial response is proportional to the transverse one, the factor between the two being roughly 
\/l -|- l/ravG^^ /G^^) : this agrees with the conjecture of |41, which is proven here to be valid 
within a fully-relativistic context. It also appears that the self-energy contribution increases in 
going from q = 500 MeV/c to 1 GeV/c, then saturates for higher values of q (thus scaling of first 
kind is fulfilled). This is due to the same particle- hole cancellation mechanism occurring in the 
electromagnetic case. 

In the longitudinal channel the self-energy gives instead a very small contribution compared 
with the MEC and vertex correlations. The effect of the SE contribution is found to increase with 
q, but always remains smaller than the one arising from the other correlations. Indeed the one-body 
longitudinal response is suppressed due to a delicate cancellation between the isoscalar and isovector 
responses [^]. Physically this occurrence reflects the fact that the electric form factor in one of 
the two vertices of diagram 18a is always very small, for both protons and neutrons. When isospin 



correlations are taken into account this balance can be disrupted, as pointed out in [51, 55, 67], 
thus yielding the large effects observed in Fig. 19. Indeed in the diagrams 18e-f a proton can be 
converted into a neutron, leading to two large couplings, Gep and Geu- Using different language, 
the MEC and VC are more effective than the SE in R\y-, since they act differently in the two 
isospin channels. Indeed the VC carries a factor -3 in the isoscalar response and +1 in the isovector 
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Figure 19: The longitudinal (left panels), transverse (central panels) and axial (right panels) PV 
responses plotted versus Dashed line: one-body contribution; dotted line: one-body+MEC+VC; 
solid line: total (including SE). 
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Figure 20: Separate contributions to the PV longitudinal (left panels), transverse (central panels) 
and axial (right panels) responses plotted versus uj. Solid: VC; dashed: seagull; dot-dashed: pion- 
in- flight. 
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Figure 21: The PV asymmetry displayed versus co for various values of the momentum transfer q 
and the scattering angle 6e- Dashed: one-body; dot-dashed: one-body-|-MEC-|-VC; solid: total. 
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one and the MEC are purely isovector, whereas the SE has almost the same impact in the two 
channels. 

This is clearly illustrated in Fig. 20, where the separate seagull (dashed), pion-in-flight (dot- 
dashed) and VC (solid) contributions are displayed. In the L channel the role of MEC is almost 
negligible, in agreement with the findings for the electromagnetic , whereas the effect of the 
vertex correlations is dominant, especially at small values of q. In the T and T' channels the balance 
between MEC and VC is similar to that occurring for the electromagnetic i?^ (see Fig. 15): the 
pion-in-flight gives the smallest contribution, particularly for large q, whereas the seagull dominates 
for all q and tends to cancel the VC contribution. Note also that the seagull and VC vanish exactly 
at the same value of w, which, for high g, coincides with the QEP. 

Since the three PV responses are not at present experimentally separable, we now explore the 
effect of the pionic physics on the asymmetry in eq. ( |154| ). In Fig. 21 we show A at various values 
of the momentum transfer q and of the electron scattering angle 9e for the free RFC (dashed) , and 
including the MEC and VC (dot-dashed) or the MEC, VC and SE (solid) contributions. Clearly the 
pionic correlations are mostly felt at low values of (left panel), where the longitudinal response 
is enhanced by the kinematical factor vl, and at low values of q, where the vertex correlations 
dominate. At high values of 9e (right panel) the asymmetry is totally insensitive to pions, because 
the effect of the SE (which gives the main contribution) cancels between the PV and PC responses 
appearing in the numerator and denominator of eq. ( |154| ). 

We thus conclude that the extraction (at large electron angles) of the axial nucleonic form factor 
Ga is almost independent of the nuclear model. On the contrary at small angles PV experiments can 
measure the strange electric content of the nucleon only if a good control of the nuclear dynamics 
is achieved, since the isospin correlations give very large effects. Conversely, interesting insight into 
the latter can in principle be gained here. Our results show that only at very large momentum 
transfer does the forward-angle asymmetry become insensitive to pionic correlations and hence 
suitable for assessing the strangeness content of the nucleon. 



4 Non-relativistic reductions 



For years most of the effects introduced by the two-body pionic currents in electron scattering 

,|70|, 



reactions have been explored assuming different types of non-relativistic reduction |28, 68 



71, U%. Not only non-relativistic wave functions have been used, but also non-relativistic current 
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operators derived from a direct Pauli reduction have been considered. Focusing on the pionic effects 
on the hadronic (e, e') response functions one has to deal with the single-nucleon electromagnetic 
(electroweak in general) current and the various two-body pionic currents discussed previously. 
Concerning the former, an improved version of the single-nucleon electromagnetic current has been 



suggested in |66|, where the expression of the current is derived as a non-relativistic expansion in 
terms of the dimensionless parameter rj = p/m, p being the three-momentum of the struck nucleon. 
In Appendix F we review this approach — which we call semi-relativistic (SR) — and compare it 
with the traditional non-relativistic reduction, where the non-relativistic expansion is performed 



with the additional assumption k = q/2m << 1 and A = uj/2m << 1. As shown in [20, ^ and 
in Appendix F, the expansion of the current to first order in the variable r/ yields quite simple 
expressions; moreover the various pieces of the relativized current differ from the traditional non- 
relativistic expressions only by multiplicative (g, a;)-dependent factors, and therefore are easy to 
implement in already existing non-relativistic models. 

The semi-relativistic form of the OB electromagnetic current operator was first checked in |66], 
where the inclusive longitudinal and transverse responses of a non-relativistic Fermi gas were found 
to agree with the exact relativistic result within a few percent if one uses relativistic kinematics 
when computing the energy of the ejected nucleon. Recently the same expansion has been tested 
with great success by comparing with the relativistic exclusive polarized responses for the ■^H(e, e'p) 



reaction at high momentum transfers [73|. This relativized current has also been applied to the 
calculation of inclusive and exclusive responses that arise in the scattering of polarized electrons 



from unpolarized |74| and polarized nuclei |7^, 77, Finally, it also has been compared with 
a fully relativistic DWIA calculation of (e, e'p) observables for |Qp = 0.8(GeV/c)2 in ||7|, A 



systematic analysis of the semi-relativistic approximation in the case of (e, e'A^) reactions has been 
presented in ||8l|| . 

Alternative expansions, in powers of the initial nucleon momentum, of the structure functions 
of nuclei have recently been proposed [15] and "recipes" to obtain a relativistic structure function 
from its non-relativistic analog by changing the scaling variable and performing an energy shift 
have been suggested. The so called three-dimensional reduction model, which includes final-state 
interactions, has been tested in the case of a deuteron target, but not for A > 3. 

The necessity of a semi-relativistic form for the current even for moderate momentum transfer 
is demonstrated in Fig. 22, where we compare the traditional non-relativistic results for the elec- 
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Figure 22: The one-body longitudinal (left panels) and transverse (right panels) responses displayed 
versus lo for various values of the momentum transfer q and of the Fermi momentum kp- Dashed: 
non-relativistic; solid: relativistic. 
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tromagnetic responses with the fully-relativistic calculation for a non-interacting system. It clearly 
appears that for low densities and momentum transfers the two approaches are equivalent, but 
that the two curves deviate from each other as q and kp increase. One of the effects of relativity is 



the shrinking of the response region [25| and is already significant at g = 400 MeV/c. This effect, 
which arises from the relativistic kinematics in the energy-conserving delta-function appearing in 
the responses, can be accounted for approximately by the replacement 

A ^ A(l + A) . (172) 

Another effect, stemming from the non-relativistic reduction of the currents, relates to the enhance- 
ment of the longitudinal response and to the reduction of the transverse one due to relativity. Such 
an effect can be mimicked by the kinematical factors k?/t (in the L channel) and t / k? (in the T 
channel), which naturally emerge from the r/ expansion illustrated in Appendix F. When included 



in the non-relativistic responses these factors, together with the prescription of eq. (172), allow one 
to reproduce the fully-relativistic responses even for very high g- values (see, for example, []66|| ). 
In this section we explore the impact of relativity on the meson-exchange currents. 

4.1 Pion exchange currents 

We first compare the fully-relativistic transverse MEC responses with the traditional non-relativistic 
calculation developed in []82|| , where the seagull p-h matrix element is evaluated analytically, while 
the pion-in-flight contribution is reduced to a one-dimensional integral. For this comparison the 
value Ftt = 1 for the ttNN form factor and the static pion propagator have been used in the 
relativistic calculation. The effect of static versus dynamic pion propagator will be discussed later 
on. 

From Fig. 23 it emerges that the two calculations give the same results for small density and 
momentum transfer. As q and kp increase we see that, apart from the difference stemming from 
the relativistic kinematics, which shrinks the response domain, the relativistic responses are smaller 
than the non-relativistic ones: this reduction amounts to about 30% for g=500 MeV/c and kF=250 
MeV/c, indicating that relativity plays an important role even for not so high g- values. 

The same curves are displayed for /ci?=250 MeV/c and g=500, 600 and 700 MeV/c in Fig. 24, 
where it is shown that the effect of relativity clearly grows with the momentum transfer. 

In Fig. 25 the relativistic MEC-correlated transverse response (dotted) is compared with the 
corresponding non-relativistic one (dot-dashed) as well as with the relativistic (solid) and non- 
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relativistic (dashed) one-body response for three values of q. The figure shows that for low values 
of q (500 MeV/c) the effects of MEC and relativity are roughly of the same size, the former acting 
mainly to the left of the QEP, the latter to the right. As q increases, the effect of relativity becomes 
dominant, pointing to the necessity of a relativistic treatment for momentum transfers larger than 
500 MeV/c. 

Finally the impact on the responses of the relativistic propagator IS.t^{K) = [K"^ — m^)~^ as 
compared with the static one aI""'"' ■*(]<;) = — (k^+m^)^"'^, which is commonly used in non-relativistic 
calculations, is explored. In Fig. 26 the pion-in- flight, seagull and total MEC contributions to 
are evaluated for (7=0.5 and 2 GeV/c using the two versions of the propagator. It appears that the 
dynamical propagator affects the pion-in-flight contribution more than the seagull term (it increases 
the latter by more than a factor 2 at q=2 GeV/c); however, the two effects tend to cancel, so that 
their net effect is not very significant. 

4.1.1 The rjp expansion 

In view of the relevance of relativistic effects illustrated above and following the ideas and methods 
developed in the case of the single-nucleon electromagnetic current operator and its non-relativistic 
reduction |66], a new semi-relativistic reduction of the MEC has been developed in pO], where 



the transferred energy and momentum have been left unexpanded while expanding only the initial 
nucleon momentum. The expressions thus obtained retain important aspects of relativity not 
included in the traditional non-relativistic MEC used throughout the literature. Here we summarize 
the basic results. 

We are interested in the evaluation of the particle-hole matrix elements {pk\j^i ^fj\kh) and their 
new semi-relativistic expressions. The resulting expansion for the MEC should be used together 



with the single-nucleon current, developed to first order in rj (see |20| and Appendix F), to set up 
the various responses. Therefore, in order to be consistent, the expansion of the MEC should also 
be performed to first order in the corresponding small quantities {r/^ = k/m,r]fi = h/m}, whereas 
{r]p = p/m,K = q/2m} are treated exactly. 

After some algebra |2^ the following semi-relativistic (SR) expressions of the MEC currents 
(referred to as NRl approximation in [pO|) are obtained: 
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Figure 23: The seagull (left panels) and pion-in-flight (right panels) contributions to the trans- 
verse response displayed versus to for various values of the momentum transfer q and of the Fermi 
momentum kp- Dashed: non-relativistic; solid: relativistic. 
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Figure 24: The seagull (left panels), pion- in-flight (central panels) and total MEC (right panels) 
contributions to the transverse response displayed versus uj for ^^^=250 MeV/c and various values 
of the momentum transfer g. Dashed: non-relativistic; solid: relativistic. 
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Figure 25: The transverse response displayed versus u for A;i?=250 MeV/c and various values 
of the momentum transfer q. Dashed: one-body non-relativistic; solid: RFG; dot-dashed: one- 

body+MEC non-relativistic; dotted: RFG+MEC relativistic. 
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Figure 26: MEC contribution to (in GeV^'^) versus uj with dynamic (solid curves) and static 
(dashed curves) pion propagator at (a) q= 0.5 and (b) 2 GeV/c. The separate pion-in-flight and 
seagull contributions are displayed. 

Seagull Current Operator. 



J°(p,k,k,p)sm 



I" 



--X's 



cr • [2ft + + (1 + T)r/fc] Xs^xl^o- ■ {rj^ - rj^^) 



(K - HY - ml 



X. 



Sh 



(173) 



2cr • K 1 



X X 



where the factor 



2(1 + T 
(1 + r) 



+ cr • (^h - Vk) - • 



X. 



Sk 



{K - Hf 



mt 



-Xs^xLfr ■ {Vk - Vh)) Xsh , 



pm 



i^Zab {tp\'^a\tk){tk\Th\th)Fl 



(174) 



(175) 



has been introduced. Note that if the terms rfj^ — {1 + T)r}i^ and 77^^ -|- (1 -|- t)t7^ are neglected 
(this approximation will be referred to as SR2) the expression for the time component is similar to 
the one obtained in the traditional non-relativistic reduction [^6| except for the factor 1/vT+T, 
which accordingly incorporates important aspects of relativity not considered in the traditional 
non-relativistic reduction. Analogously, in the space component, if the terms {k ■ r7^)/[2(l -|- r)] 
and T{(T-r]j^) are neglected {SR2 approximation), the traditional non-relativistic expression is 



simply recovered, except for the factors l/y/l + r and ^/T^\^ that multiply the contributions given 
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Figure 27: Seagull current matrix element -fC^ - see eq. ( |176| ) - for q=\ GeV/c and kp= 250 
MeV/c. The kinematics for the hole are h=175 MeV/c and <j)h = 0. First column: spin (1/2, 1/2) 
component; second column: spin (1/2,-1/2) component. Solid: fully-relativistic; dashed: SRI 
approximation; dot-dashed: traditional non-relativistic; dotted: SR2 approximation. 
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Figure 28: Pion-in flight current matrix element K'j^ for q=l GeV/c and /ci?= 250 MeV/c. The 
kinematics for the hole are /i=175 MeV/c and 0/i = 0. First column: spin (1/2, 1/2) component; 
second column: spin (1/2,-1/2) component. Solid: fully-relativistic; dashed: semi-relativistic 
approximation; dot-dashed: traditional non-relativistic. 
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by the two diagrams involved. Thus, as in the case of the time component, also here important 
relativistic effects are simply accounted for by these multiplicative factors. 

To illustrate this point we plot in Fig. 27 the current matrix element i^^ defined through 

f E £i;(p.'<.>'.i.) = --Ji=jf K;(,,.,h)i:^ (176) 

k<kp ^ yi^/p-fch-C/k 

for q=l GeV/c, /i=175 MeV/c and (j)h=0. The curves represent the fully-relativistic result (solid), 
the traditional non-relativistic approximation, including relativistic kinematics through eq. ( |172| ) 
(dot-dashed), the SRI approximation of eqs. ( [173 , 174|) (dashed) and the SR2 approximation (dot- 



ted). Only the relevant components are shown, the other vanishing for symmetry reasons (see |2^ 
for details). It clearly appears that, while the traditional non-relativistic reduction, although cor- 
rected by the replacement A A(l -|- A), fails to reproduce the exact results by roughly 10 to 
20% (this deviation increasing with q, as shown in ||2^), both the SRI and SR2 approaches yield 
excellent agreement with the fully-relativistic current. 

Finally, we examine the limit rjp ^ 0, since this provides some understanding of how the MEG 
effects are expected to evolve in going from light {r]F very small) to heavy nuclei {rfp = 0.29). In 
this limit the seagull current simply reduces to 

hm jO(p,k,k,h) =0 . (177) 

This is because the time component of the seagull current is of first order in the small momenta 
involved or, equivalently, it is O^-qp). On the contrary, the vector component in the limit ?7i? — > 
becomes 



lim i',(p,k,k,h) = ^^ xL ^"' J^2^"'\'''^ 



t .. 

^^-^2 > (178) 



which shows that the space components of the seagull current are 0(1) and contribute even for 
nucleons at rest, as happens for the charge and magnetization pieces of the one-body current. 

Pion-in-flight Current Operator. 

Keeping only linear terms in the small momenta, one obtains for the semi-relativistic pion-in- 
flight current 

j,(p,k,k,h),« = -^^^m xl^ _ _ _ _ ^^^^ xs^n . (180) 
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Again, these expressions are similar to the traditional non-relativistic currents |24] except for the 
common factor \l \J\ + r, which includes important aspects of relativity not taken into account in 
the traditional non-relativistic reduction. Note that the space component of the pionic current is, 
in leading order, purely longitudinal; its transverse components are in fact of second order in r]p. 

In Fig. 28 we display the current matrix element -fC^, defined analogously to eq. ( |176D , for q=\ 
GeV/c, h=175 MeV/c and 0/1=0. As for the seagull, the fully-relativistic result (solid) is very well 
approximated by the SR prescription (dashed), whereas the traditional non-relativistic approach 
with relativistic kinematics (dot-dashed) deviates from the exact result by 10-20%. 

Finally the limit rjp ^ implies that j^(p,k, k, h) = 0, since all components of the pionic 
current are 0{r]p) in the expansion. 

Summarizing, the r] expansion shows that relativistic kinematics can be very easily implemented 
in MEC semi-relativistic calculations by applying the prescription A — > A(l + A) and by multiplying 
the exchange currents by the kinematical factors indicated in eqs. (|173| - |180|) . 



4.1.2 Results for the responses using the relativized MEC model 

In this section we discuss the validity of the relativizing prescriptions introduced above, when they 
are implemented in a traditional non-relativistic model of the reaction. We begin with the non- 
relativistic Fermi gas of [^2|, which includes also MEC in the transverse response. One of the 
advantages of this model is that the integral over the Fermi sea appearing in the seagull matrix 
elements can be performed analytically, while the pion-in-flight is reduced to an one-dimensional 
integral. 

Next we will relativize this model by implementing relativistic kinematics through the substi- 
tution A A(l -|- A) in all places except in the nucleon and pion form factors F{q, w), which should 
be evaluated at the correct cj-value. Second, we use the new semi-relativistic expansion of the 
electromagnetic OB-I-MEC operators in powers of r]. For the OB operators we use the following 
expressions (see Appendix F): 

POB = -^Ge+1 J== {KXr])-CT (181) 

V'T Vl + r 

Jos = — [iGMio- X k) + GeVt] ■ (182) 

K 

Note that near the QEP it makes little difference to use the factors 1 -|- r or k^/t. In these 
factors lies the main difference with the traditional non-relativistic charge and transverse current 
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operators. Note that in addition we include a first-order spin-orbit term in the charge operator. 
The contribution of this term is small in the longitudinal unpolarized response, since its interference 
with the leading order is exactly zero in PWIA and hence it gives a negligible contribution of second 
order in rj. However one should be careful in including this term in the more complex cases in which 
there is an interference TL response (when the nucleus is polarized or in the exclusive reactions 
{e,e'p), see |7^, Q), and where it gives a significant contribution, since in this response the 
leading order is zero in PWIA. 

In the case of the MEC we use the following simplified prescription to relativize transverse 
operators: 

JIiEC = / I , JMEC,nonrel ) (183) 

namely we introduce a factor 1/^1 + t to take into account relativistic corrections coming from 
the free Dirac spinors. Note that in the case of the seagull we have neglected a further correction 
factor 1 -|-r in the hole part of the seagull current. However here we choose to present results with 
the above simplified version of the transverse current, since it is the easier to implement in already 
existing models of the reaction; otherwise one has to identify the different pieces of the operator, 
which may be difficult. Furthermore, this correction is not of much importance, its main effect 
being to correct slightly the position of the zero in the seagull response. Be it as it may this ad hoc 
prescription for the seagull current is supported by the quality of the results shown below. 



Results for the *"Ca nucleus for q = 500 MeV/c and q = 1000 MeV/c are shown in Figs. [29| 



and 30, respectively. In the upper part of these figures we show the one-body (OB) separated 
longitudinal and transverse responses. The solid lines are the exact relativistic results. These are 
very different from the traditional non-relativistic results shown with dashed lines. Note that the 
same nucleon form factors and the same kp = 237 MeV/c are used in both calculations. If we 
include relativistic kinematics, then we obtain the dotted lines, which are still different from the 
exact result, even if now the region where the response is nonzero is similar to the relativistic case. 
Finally, using in addition the new semi-relativistic corrections (factors k/^/t in the charge and 
\/t/k in the current) we obtain the relativistic approximation shown with dot-dashed lines, which 
is very similar to the exact result. Hence we can safely say that the new expansion of the OB 
current is very good, giving essentially the exact answer. 



The case of the MEC transverse responses is shown in the lower part of Figs. M and 30. There 



we show the separate contribution of seagull and pion-in-flight currents to the transverse response 
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Figure 29: One-body longitudinal and transverse response functions (top panels), and transverse responses 
of interference between MEG and OB currents (bottom panels) for q = 500 MeV/c and kp = 237 MeV/c. 
Solid: exact relativistic results with static propagator and without ttN form factor. The rest of the curves 
have been computed using the non-relativistic Fermi gas model, with or without relativistic corrections. 
Dashed: traditional non-relativistic results. Dotted: including relativistic kinematics in the non-relativistic 
calculations. Dot-dashed: including in addition the new expansion of the OB currents. Double-dashed: 
including in addition a correction to the MEG operators with a factor 1 / y/1 + r. 

(interference with the OB current). Again we show with solid Unes the exact relativistic results, 
and with dashed lines the traditional non-relativistic results. If again we include the relativistic 
kinematics we obtain the dotted lines. With dot-dashed lines wc display results which include in 
addition the relativistic correction to the OB current, amounting to a factor l/vT+T ~ y^/n. 
This correction produces a small reduction of the responses. Finally, with double-dashed lines we 
show the results computed using in addition the relativistic corrections in the MEC, which amounts 
to another factor l/\/r+T. This correction produces a further reduction of the responses, giving 
a result which is closer to the exact one. 

From these results it appears that our expansion of MEC currents is not as good (at least 
fractionally) as the OB expansion. This is likely related to the fact that the OB currents have been 
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Figure 30: The same as Fig. |9[ now for q ^ 1000 MeV/c. 

expanded in powers of r/ = h/M, where h is tlie momentum of the hole, and they are exact, by 
construction, for h = 0. However in the case of the MEC there is another variable in the expansion: 
the momentum of the second hole h' /M, which is small, but is never zero and is being integrated 
up to kp. Therefore our currents are not constructed to agree with the exact ones for h = 0. 

On the other hand, we have explicitly showed before that in the limit q ^ and kp — > 0, the 
relativistic and non-relativistic results agree. This is also the case for the present results of the 
relativized currents, as it is illustrated in Fig. There we show the seagull and pionic responses 
for several small values oi q = 100, . . . , 500 MeV/c and for kp = q/2. With solid lines we show the 
exact relativistic results, while with dashed lines we show the traditional non-relativistic results. 
Finally we also show with dotted lines the results using the present semi-relativized approach. It 
is seen that the last are always much closer to the exact result than the non-relativistic ones, and 
that they converge faster to the exact results. 

Better agreement between the exact and the relativized models for the MEC responses is also 
expected in the limit rjp ^ in the quasielastic peak, since in this case both momenta h and h' are 
forced to be small, which are the conditions assumed in our expansion. Results for the transverse 
MEC responses in this limit are shown in Fig. ^ for q = 1000 MeV/c and for two values oi kp = 50 
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Seagull 



Pion in flight 




Figure 31: Seagull and pionic responses computed for several values of q and kp- Solid lines: exact 
relativistic results. Dashed lines: non-relativistic results. Dotted lines: approximated semi-relativistic results 
using relativistic kinematics and relativized currents. Static propagators without a ttNN form factor have 
been used here. 
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Figure 32: MEC-OB transverse responses for for q = 1000 MeV/c and kp = 50 MeV/c. Solid: exact 
relativistic results with static propagator and without a ttN form factor. Dashed: non-relativistic results, 
but including relativistic kinematics. Dotted: using in addition the relativized currents. 

and 25 MeV/c. With solid lines we show the exact relativistic results, while with dashed lines we 
display the non-relativistic ones, but including relativistic kinematics. Finally, the dotted lines 
correspond to the semi-relativized results, which fully agree with the exact ones in the case of the 
seagull current, while in the pionic case there is still a difference between the two calculations. 
However this is not very important in this limit, since the transverse pionic contribution which we 
are computing here is of second order in the hole momenta, and so is negligible compared with the 
seagull one, as can be seen in the figure. 

In order to improve the present results for the MEC responses one should look for an expansion 
of the MEC in the form 

Jip, h) = ^^]f'^^h ^,^UP, h) , (184) 
where iV(g, a;, kp) is an appropriate normalization factor defined by requiring 

lim-^ = l, (185) 
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i.e., the coincidence between the relativistic and the approximate results at the quasielastic peak. 
Obviously the factor N{q,uj, kp) is a function of kp also, since an integral over the Fermi sphere is 
implicit in the definition of the MEG in the Ip-lh channel, and it can be written in the form 

N{q, u;,kp) = TTTT /-'^'^;"^ . (186) 

•""non rel W' ^) 

A simple approximation for this function is not easy to obtain, since it requires the knowledge of 
the exact relativistic answer. 

Despite these difficulties, the quality of the OB expansion plus the approximated improvement of 
the MEG currents obtained in the present expansion are good enough to ensure a quite satisfactory 
description of the exact relativistic transverse response using the relativized OB plus MEG operators 



altogether with relativistic kinematics. This is shown in Fig. |3^, where we show the total transverse 
response, including OB+MEG operators, for q = 500, 1000, 2000, and 3000 MeV/c. The sohd lines 
are the exact relativistic result. Again with dashed lines we display the traditional non-relativistic 
results, which together with relativistic kinematics give the dotted lines. Finally, with dot-dashed 
lines we show the results using the semi-relativized OB+MEG currents. The agreement between 
the two models is quite good even for very high q, since the major part of the relativistic effects 
is included in the approximated model, and therefore these currents are very appropriate and easy 
to implement in already existing non-relativistic models of the reaction. 

4.1.3 Comparison with the traditional relativistic corrections 

Here we discuss the reasons why the present expansion of electromagnetic operators is preferable 
to other kinds of relativistic corrections existing in the literature. The most common of these is 



the Darwin- Foldy correction to the charge operator of the nucleus |0, q^, 84, 85]. This correction 
is usually derived from a Foldy- Wouthuysen transformation but can also be obtained from a 
Pauli reduction of the spin matrix element (we do not write the spin indices for simplicity) 

J^(p, h) = UE{p)T^'{Q)uE{h) = ue{p) (^Fi7^ + i^cT^-Q,^ .^^(h) , (187) 

where we use the sub-index E to denote the spinors normalized to u^p{p)uE{p) = 1, i.e., namely 
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Figure 33: Total transverse response function of ^°Ca including MEC for several values of the momentum 
transfer, and for kp = 237 MeV/c. Solid: exact relativistic results. The rest of the curves have been computed 
using the non-relativistic Fermi gas model, with or without relativistic corrections. Dashed: traditional non- 
relativistic results. Dotted: including relativistic kinematics in the non-relativistic calculations. Dot-dashed: 
including in addition the new expansion of the OB+MEC currents. The relativistic calculations include 
a dynamical propagator and TriV form factor, while the non-relativistic calculations do not include these 
corrections. 
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This is in contrast to the Bjorken and Drell spinor normaUzation used in the present work, where 
spinors are normaUzed to n^(p)u(p) = E/m. The relation between the two sets of spinors obviously 
is 



«(P) = \ —ue{p) ■ 
m 



(189) 



Of course both formalisms based on the different spinor sets ue{v) or ^(p) should give the same 
results for the observable quantities. For instance, if we like to think in terms of wave functions, 
for a nucleon in a box of volume V this would be 

1 



m 



UE(p)e-f-. 



(190) 



This means that observables (expectation values, probabilities or cross sections) computed using 
the Bjorken and Drell normalization always contain additional phase-space factors m/E, while 
these factors do not appear explicitly if one uses the i?-scheme, since they are already included 
inside the spinors UEip)- 

As an example let us consider the case of the longitudinal response function for protons 

- EE ^(^p-^h-a;)|(p/i-^|p(q)|F)P 
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ph SpSfi 
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P(P, h)V(P,h) 



d^h S{Ep -Ei,-u>) 
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-Tr 



p(P,h)V(P,h) 



(191) 



where in the last line p = h + q, and we have used the replacement J2h — *■ / with 
V/{2Tr)^ = SZ/Sirkp. Note that we use the Bjorken and Drell normalization and so the energy 
denominators appear explicitly. The charge matrix element used here is the fully-relativistic one 

p(p,h) = «(p)rO(g)u(h). 

The interesting (and crucial) point is that the energy denominator Ep cancels out when we 
perform the integral over cos 6 — the angle between h and q — using the energy-conserving delta 
function. In fact, from 



E' 



p2 + ^2 ^ j^2 _^ q2 _^ cos e + rri^ 



(192) 



we have EpdEp 



E, 



hqd cos 9. Therefore the angle 9 becomes fixed by the energy conservation 
p — Eh + to and we obtain 

3Z r^p 



8nkl 



F JO 
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p(p,h)V(P,h) 



(193) 



This expression has to be compared with the the non-relativistic response function, which can 
be computed by repeating the steps above using instead the non-relativistic energies Cp = p^/2m, 
i.e., 

( 3Z 



non rel 



h<kp 



fi'/i 5(ep - eh - c^)TV p(p,h)V(p,h) 



(194) 



where again p = h + q, no energy denominators appear and the non-relativistic charge operator is 
used. 

Now the integral over cos 9 can again be performed using the non-relativistic identity 
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2m 



2 q2 + 2hq cos / 



2m 



(195) 



Hence de-p = ^ cos 9 and the integral over the new variable ep can be performed. The latter is 



eh + w. We obtain 



3Z m r^F 



non rcl 



hdh 



2tv 



#Tr p(p,h)V(p,h) 



(196) 



STTk\ q Jo 

which has formally the same structure as eq. (|193| ) with the exception of the factor m/E\^ ~ 1 
included there. Therefore, the relativistic response can be reproduced using a non-relativistic model 
if we introduce in the non-relativistic response in eq. (|196| ) a good approximation for /9(p,h), and 
in addition we use relativistic kinematics, i.e., we use the relativistic relation between cos 9 and w. 
This can be approximately accomplished starting from eq. (|192|). Indeed we have 



h^ + q^ + 2hq cos i 



(^h + w 



- = h^ +u;^ + 2Eh + co 



~ h"^ +uj'^ + 2mijj + 2e\^oj 



2m \ijO + e^ + 



ijj 



2eh^ \ 



2m J 



2m \ uj + eh + 



UJ 



2 \ 



(197) 



2m I ' 

where we have neglected the term ehuj/m = 0{h^ /m^). Comparing with the non-relativistic relation 
( |195| ) we see that the relativistic one can be approximately obtained with the replacement u) — > w(l-|- 
uj/2m). The validity of this approximate method of relativization was demonstrated numerically 
in the last sections. 
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If instead we use the other spinor normahzation ue{p) in the relativistic model, then the 
matrix element is /9^(p, h) and now there are no explicit energy denominators in the expression of 
the response, which hence becomes, after integration over 9, 

RL{q^u^) = —^- j hdh I d(l)-^TT\pEip,h)^PE{p,h)\ ■ (198) 



F q Jo Jo 



m 



Comparing with eq. ( 196 ) we see that if we start with a non-relativistic model and use an approxi- 
mate non-relativistic form for peIp, h), as happens with the Darwin-Foldy correction, an additional 
factor Ep/m is needed in order to reproduce the relativistic response. For this reason, a careless 
introduction of relativistic corrections alone in non-relativistic models can produce incorrect results. 

4.2 Pionic correlations 

The analysis carried out in the previous section for the MEC could in principle be extended to the 
correlation current. However in this case the calculation becomes extremely cumbersome and has 
not yet been performed. A semi-relativistic calculation of the vertex and self-energy correlations 



has been carried out in |M, 55, p6[, where the relativistic energy-conserving delta function has 



been accounted for via the replacement in eq. ( 172 ) and the form factors in the two-body current 



have been modified to implement relativistic effects. The response functions so obtained are in 
qualitative agreement with the fully-relativistic ones for not too high q. However, for high values 
of g a careful treatment of the relativistic effects is needed. 

In what follows we briefly examine the non-relativistic limit of the vertex correlations and self- 
energy diagrams in order to bring to light some differences with respect to the fully-relativistic 
case. 



The non-relativistic leading order of the pionic correlation currents in eqs. (p7|,28) is obtained 
by using the following prescriptions, valid in the static limit: 



E^k =^ m (199) 

(200) 
(201) 



75 ^ ~ cr • k (200) 
1 1 



2 



Sf{P) SnriP) = ^—r ■ (202) 

The electromagnetic form factor r'^(Q) is also replaced by T!^^{q), representing the usual non- 



relativistic one-body current acting over bi-spinors |2C, pj]. Using the above relations and perform- 



ing the sums over spin and isospin indices, the VC and SE current matrix elements read 
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J^ciPMnr = E { (ifj +^2 Snr{K + g)r.C.(Q)r,a • (k - h) 

+ ^.(p-kKr^:,(Q)r,5„,(i^-Q)-^-[H_iiI^\^ (203) 

(p - k)"^ + J 

and 

j^g(p, h)„, xl, [^nr{v)Snr{P)Kr{Q) + r^.(Q)5„r(i^)S„.(h)] Xs, , (204) 

where Xsp and Xs^ are two-components spinors. The non-relativistic self-energy function is given 
by 

(p) = (pVT„; = ^■""■'i' ■ 

The above expressions coincide with the traditional non-relativistic currents used in the literature. 

With the self-energy in eq. ( p05[ ) one can then construct the non-relativistic Fock nucleon 
propagator 

5f/(po,p) = ^ + ^S„,(p) + • • • = . (206) 

As is well-known, this is a meromorphic function whose simple pole again defines the new energy 
of the nucleon in the medium, namely 

enr(p) = ^ + S„,,(p) , (207) 

since S„r(p) is a function only of p. 

Since the non-relativistic self-energy function in eq. ( 205D does not depend on spin, the nucleon 



wave functions are not modified in the medium. In fact the corresponding Schrodinger equation in 

momentum space, including the self-energy, is simply given by 
" ^2 



^ + S„r(p) 
2m 



4>nrip) =PO</'nr(p) , (208) 



with the bi-spinor (pnrip) corresponding to the eigenvalue po = e„r(p)- 

The non-relativistic analysis of the nucleon self-energy current is much simpler than its 
relativistic counterpart. Indeed, in the former the self-consistency is immediately achieved because 
the nucleon wave functions are not modified by the self-energy interaction and thus the first iteration 
of the "Hartree-Fock" equations already provides the exact energy. By contrast, in the relativistic 
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framework the spin dependence of the self-energy [^7|| modifies the Dirac-spinors, inducing an 
enhancement of the lower components. Moreover, the field-strength renormalization constant, 
namely the residue of the nucleon propagator in eq. ( ^061) at the pole, in the non-relativistic case 
is just unity. Hence the enhancement of the lower components and the spinors' field strength 
renormalization are genuine relativistic effects absent in a non-relativistic analysis where only the 
energy-momentum relation in the medium is altered by the self-energy diagrams. We have shown in 
Section that the two above-mentioned relativistic signatures can be incorporated as new pieces 
in the electromagnetic current acting over free spinors. 

5 Conclusions 

Our goal in these studies has been to explore some of the ingredients that enter at high energies 
where relativistic effects become important in attempting to model the nuclear response functions 
for inclusive quasielastic electron scattering. The full problem of accounting for relativistic dynamics 
in nuclear physics is a daunting one and far from being solved While in many papers 

it appears that a reasonable level of understanding has been reached |||, since the 

basic trends seen in the data are reproduced, closer scrutiny reveals a different situation. It is 
not only that contributions left out in various analyses are far from being small, but, even more 
serious, fundamental physics principles (Lorentz covariance, gauge invariance and unitarity) turn 
out patently to be violated. Thus the successes in reproducing the experiments often reflect more 
an adjusting of parameters than a real understanding of the physics involved in the quasielastic 
regime. 

Our approach has been less to use a highly elaborated non-relativistic model whose failings 
are expected at the outset than to employ a simple model in which the important consequences 
of relativity are hopefully present. For this we have begun with the relativistic Fermi gas as our 
starting point |9^. This approach is motivated by several critical features of the model, namely, 
that it is Lorentz covariant, that it allows the implementation of gauge invariance and that is it 
simple enough to be tractable and yet not obviously lacking at least for the quasielastic responses 
for which it is designed. Clearly it is not an appropriate way to proceed if near-Fermi-surface 
physics is the goal and this regime is not our focus. 

With these as basic motivations in a series of papers we have explored the consequences of 
having a Lorentz covariant model. In particular in pO] we attempted to approximate the full 
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theory by identifying a dimensionless variable that is smah enough to be suitable in setting up 
a semi-relativistic expansion of the responses (namely the momentum of a nucleon lying below 



the Fermi surface compared with its mass). In contrast, in very recent work |1£] no expansion 
whatsoever is involved and the theory is now fully-relativistic. 

Our treatment proceeds in terms of nucleonic and mesonic degrees of freedom (the latter viewed 
both as force and current carriers). As our aim is to study the quasielastic regime where the longest- 
range hadronic ingredients may be expected to be dominant, we focus on pions; studies using a 
larger set of hadrons can be undertaken and some steps have already been taken by us in that 
direction. In our model, the pions are dealt with to first order in a perturbative framework, since 
their effects on the free responses of the RFG are not expected to be too disruptive. 

Gauge invariance is a fundamental property we have also addressed in very recent work [^]. We 
now understand how the continuity equation is satisfied order by order in perturbation theory. We 
have succeeded in showing that the continuity equation for the one-body (single-nucleon) and the 
two-body (MEG and correlations) currents is fulfilled, implying that our approach deals consistently 
with both forces and currents. 

Given the point in our understanding of the quasielastic responses at relatively high energies, 
we have been motivated to provide a comprehensive set of discussions of progress made so far. In 



particular, as a more in-depth presentation of the analysis carried out in |19], where we first stud- 
ied the fully-relativistic set of one-pion-exchange operators that contribute to the electromagnetic 
responses of nuclei in the Ip-lh channel, in the present work we have gone further to answer the 
question of whether or not a finite OPE self-energy current in nuclear matter exists. Indeed we have 
proven that the latter can be obtained through a renormalization of the Ip-lh excitation vertex 
with a Fock self-energy insertion in the particle or in the hole line. In these diagrams were 
shown to diverge but, at the same time, to be crucial to preserve gauge invariance. To overcome this 
impasse in that work we abandoned the notion of current operators, using instead the polarization 
propagator for computation of these diagrams. Indeed the double pole appearing in the self-energy 
polarization propagator can be dealt with employing the derivative of the nucleon propagator. 

In assessing the role of the pions in the electromagnetic nuclear responses, the MEG are not the 
only contributions that arise in first-order perturbation theory. In fact the pionic correlations are 
intimately linked to MEG through the continuity equation and, as we have seen, only when the full 
set of Feynman diagrams with one pion-exchange is considered can one expect gauge invariance to 
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be fulfilled. 

Since all of these ingredients are required for a consistent theory, a question we have addressed 
in this paper is whether or not a unified treatment based on current operators at the level of 
the OPE can be used even for the self-energy contribution. We succeeded in achieving this goal 
introducing a new ingredient: the first-order correction to the wave function and energy of the 
nucleon in the medium, which is modified by the interaction with the other nucleons. Indeed the 
iteration of the self-energy diagrams generates a "dressed" propagator in the medium. By the 
same token the self-energy generates "dressed" or "renormalized" wave functions in the medium, 
solutions of an in-medium Dirac equation, where the self-energy plays the role of a mean relativistic 
potential. This equation also provides the dispersion relation linking the energy and momentum of 
the nucleon in the medium. Importantly, the new spinors should be multiplied by a renormalization 
function ■\/^(p). 

As the self-energy is generated by the interaction of a nucleon with the other nucleons in the 
medium, the solutions of the new Dirac equation should be used as input to re-compute the self- 
energy and so on. The exact answer is obtained through a self-consistent procedure. In this 
paper, however, we have just considered the first iteration: we have thus computed the self-energy 
current confining ourselves to first-order corrections to the energy and spinors - or, equivalently, 
to corrections linear in the self-energy - which correspond to diagrams with only one pionic line, 
in order to be consistent with the MEG and vertex correlation currents. 

Notably in the first-order expansion of the renormalized spinors two new elements with respect 
to the non-relativistic approach emerge, one arising from the negative-energy components in the 
wave function produced by the interaction, the other from the renormalization function \J Z2(p). 
These two elements can be combined in a new renormalized self-energy current, J^g^;, acting over 
free spinors, and, together with renormalized self-energies, lead to the same self-energy contribution 
of The introduction of renormalized energies produces a shift of the response function. Our 
results for the response functions for typical kinematics show that the negative energy components 
constitute a correction to the total self-energy contribution of roughly 10-20%, whereas the renor- 
malization function for OPE is small, yet necessary if gauge invariance is to be fulfilled exactly. 
Moreover, while at low momentum transfers both particle and hole contributions play a role in the 
response, at high q only the hole contribution survives. Finally, the self-energy contribution to the 
response functions is comparable in size to the one arising from the MEG and vertex correlations. 
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These formal developments have been discussed at length in the present work; and not to 
interrupt the flow of the arguments unduly some details have been placed in a series of Appendices. 

In the remainder of the article we have presented some typical results, both for parity-conserving 
and parity-violating quasielastic electron scattering. Briefly we have found the following: we have 
found that the MEC contributions are small enough to be well handled in first order. In particular 
both the pion-in-flight and seagull contributions are very small in the L channel where the virtual 
photon exchanged between the electron and the Fermi gas couples to the charge of the pion, 
implying as expected that the MEC only marginally affect the Coulomb sum rule. In contrast in 
the T channel the MEC are more significant. There the seagull contribution dominates, and one 
sees that the MEC contribution does not vanish when q increases. In the scaling behaviors of 
the MEC were also explored in detail: in summary it was seen that they break scaling of the second 
kind everywhere, but, while breaking scaling of the first kind at modest momentum transfers, tend 
to successful first-kind scaling behavior at sufficiently high values of q. 

The correlation contribution arising from the vertex corrections (VC) display a different pattern: 
the L channel dominates over the T channel by an amount of roughly 3:1. Thus the longitudinal 
response effectively picks up only these correlation contributions, since the MEC effects are so small 
there, and the former contribute to the total at roughly the 10-15% level. Indeed, were these to 
be the only contributions needed in addition to the REG response itself, then we would expect the 
total to shift in uj. Note that, since the correlation contributions are roughly symmetrical about 
the quasielastic peak, their impact on the Coulomb sum rule should be very small, perhaps only at 
the few percent level. The correlation contribution to is similar to the MEC contribution, but 
is smaller, roughly 1/2 the size of the latter; since the two are of opposite sign, they tend to cancel 
and thus the total is similar to the MEC contribution but is cut down by a roughly factor of two. 

In summary, the total contribution (the sum of Ip-lh MEC + Ip-lh correlations) to be added to 
the REG response (1) is not insignificant, (2) is Lorentz covariant /gauge invariant and interestingly 
(3) does not go away as q becomes very large. 
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A Gauge invariance of two-body currents 

In this Appendix we prove that the total two-body current is gauge invariant at the level of the two- 
body matrix elements in free space. We start by evaluating the contraction of the four-momentum 
transfer with the correlation current icor(p'i> P2' Pi' Ps)- It can be written as 

f-, 1 
QMicor(p'i> P2> Pi, P2) = — ^i(p'i)-ra75 2 + (1^ 2) (209) 

with M.a given by 

Ma = n(p2) [ra75 ^iSf{P2 + Q) fQFi + Fi (QSpiP^ - Q)Ta^5 ^1] U{P2) , (210) 

where we have used the relation QijT'^{Q) = Fi{Q) Q. After some algebra, involving the nucleon 
propagator and the Dirac spinors, M.a can be further simplified leading to 

Ma = n(P2) [ra75 ^li^i - i^ira75 f^l] U{V2) = U(P2) [Ta, ^i] 75 ^in(p2) . (211) 

To evaluate the commutator [ra,-Fi] we now decompose the nucleon form factor into its isoscalar 
and isovector pieces, Fi = ^ (^Ff + F^ts^ ■ Then 

[Ta,Fi] = -iFYe3abTb , (212) 

which entails the automatic conservation of the tt^ exchange current (a=3). Using eq. ( ^121) we can 
recast Ma as follows 

Ma = -iFYe3abu{P2)ni5 ^in{p2) . (213) 
Hence the divergence of the two-body correlation current matrix element can finally be written as 

j2 pV 
<5/iJcor(p'l'P2,Pl'P2) = 2m^Ze3afeM(p'i)ra75^i(pi^ 



Kf — 

X u(p^)rb75(^ + 2m)n(p2) + (1^2) . (214) 

The divergence of the seagull and pion-in-flight two-body current matrix elements can also be 
calculated in a straightforward way. The final result reads 

j!2 pV 

(3Mi^(p'i,P2>Pi>P2) = -2m— ie3a6n(p'i)ra75n(pi)— 2"^^ — - 

r^l — 

X n(P2)r675 ^«(P2) + (1^2) (215) 

QmJp(P1'P2,P1>P2) = 4m — lesab^Tr 



ml {Kl-ml){Kl-ml) 
X M(p'i)ra75n(pi)u(p2)Tb75u(p2) . (216) 
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Then, by summing up the contributions given by the correlation (eq. ( pl4| )) and seagull (eq. ( plSD ) 
currents and writing the four-momentum transfer as = {Ki + K2)^, we finally obtain 

Ql^ [jcor(p'l>P2>Pl>P2) +i^(p'l,P2,Pl,P2)] 

= Am^ ie^ab fj^2 2\(l-2 ^ 2^ '"(Pl)'^«^5^(PlMP2)^b75^(P2) , (217) 



which cancels exactly the contribution of pion-in-flight current in eq. (216) provided the electro- 
magnetic pion form factor is chosen to be i^V = . 

B Gauge invariance of the two-body current p-h matrix elements 

Following the study of gauge invariance at the level of the free-space particle-particle matrix ele- 
ments, here we extend the analysis to the particle-hole channel, deriving the contribution to the 
continuity equation of the isoscalar and isovector SE, VC and MEC particle- hole matrix elements. 
We start by evaluating the divergence of the correlation particle-hole matrix element ^^^^.(p, h) for 
the SE and VC contributions; next we address the MEC p-h matrix elements. 

• Self energy (SE) 
From eqs. (|3l| , ^2|) we get 

n -u 3/^ ^ m (/^- ^)(^-m)(f- ^) 

n ^ 3/^ ^ m ,^, (/C- ff)(^-m)(^- fi) 

Note that Fi cannot be taken out of the matrix element, since it acts on the isospinors. Now from 
the relations 

Sf{P) ^u{h) = n(h) (220) 

u{jy)^SF{H) = -u(p) (221) 

7l(p)(/'- ^)(^-m) = 2mu{v){P- f^) (222) 

{f[ -m){^- P)u{\i) = -2m(^-m)M(h) (223) 
the following expressions are derived: 

k^F {P-KY - mi 

n -u 3/2 m (^- P){^-m) 
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• Vertex correlations (VC) 

From eqs. (HQ the four-diver ffence of the VC matrix element is found to be 



Q-B: 



J2 TT^PHi^- P)SF{K + Q)TaFiTa 



^ -m 



'{K-Hf-ml 



n(h) 



E 



m 



fC — m 



E^^'^'{P-Ky-ml 



(226) 
(227) 



We now exploit the identities 

Sf{K + Q) ^{^ + m) = +{^ + m) 



to get finally 



^ ^ f ^ m ^ ^ ^ (K-H)(K-m) , 

Q • ^ = ^^(pKi-i^ /g,_%^^_^.^ n(h) 

n « \- "^-r ^ P (^-m)(f- ^) 

If the expressions ( p24 , p25 ,23C, 231| ) are split into their isoscalar and isovector parts, 
in Section 2.2.2 , we get 



^ K<kp «v \ / IT 

3/' fV ^ li)T3 ... 



P Vm 



Q ■ 



n pv- V- , (^- /f)(^-m)T3 

^< ^k,^^ {K-HY-ml ^^^^ 



= +■ 



= + 



V^/r E^ ^^(P) '"iKrH f-ml + -3..r.r.).(h) 
3/2 



+ 



3/ p5 



From these relations we learn that: 



f j,V V- , (^-m)(f- ^) , 

E ^«(P)^^3^)237;^(-3 + ^^3a.r.r,)n(h) . 



(228) 
(229) 

(230) 
(231) 
as illustrated 

(232) 
(233) 
(234) 
(235) 
(236) 
(237) 
(238) 
(239) 
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• In the isoscalar channel the self-energy and vertex contributions cancel 



Q-nf^ + Q- B^^^ = Q.n!~p + Q- = 0. 



(240) 



This differs from the non-relativistic result |24|, where the self-energy is by itself gauge in- 
variant. 



• In the isovector channel we get 



Q 



^1 ^^3a.E^^^(p) ^(h) 



2/2 



(241) 



y^2^i ^^Sab^l. ;^^i(p) Ti? TJ^1 ^1 ^(h)- (242) 



k<fci 



{K - Hy - ml 



These expressions, using the Dirac equations /fn(h) = mn(h) and u{p) f = mu{p), can be further 
simplified to yield the following four-divergence of the correlation current 



Q ■ jcoriP, h) = . [wf ) + e(^) + + .F(^) 



2/2 



k<kF 



K -P-m^ 



K-H -m'^ 



Vml ^^^'^ - |(p _ K)2 _ ^2 {K-HY-ml\ ^''''^^^ ' 



(243) 



This contribution is exactly canceled by that of the MEC (seagull and pion-in-flight) as we illustrate 
in what follows. 



MEC 



Using the expressions given in in eqs. (25,26) for the p-h matrix elements corresponding to the 
seagull and pion-in-flight currents the associated four-divergences are found to be 



Q • js(p,h) 



k.<kp 



(^- m) 



+ 



{P-KY-ml [K-Hf-ml 



n(h) (244) 



_ 2mf 
Vml 



m {Q'^ + 2H -Q -2K -Q) 

[{P - KY - ml] [{K - HY - ml 



:u{p)Ta{W -m)nu{^). (245) 



Exploiting the Dirac equation and after some algebra the above can be recast as follows 



Q ■ jMEciP, h) = Q • js(p, h) + Q • jp(p, h) 

m 



2/2 ^. m J K-P-m^ 



K-H -m^ 



Tbu(h) (246) 



We have thus proven that the correlation and MEC p-h matrix elements satisfy current conservation, 
i-e., Q ■ jcor(p,h) + Q ■ jMEciP,ii) = 0. 
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C Polarization propagator with nucleon self-energy 



Here we evaluate the Feynman diagrams for the polarization propagator with self-energy (S) inser- 
tions in the particle and hole lines, depicted in Fig. 6, subdiagrams (g) and (h). From the general 
Feynman rules for the polarization propagator [^] we have 



-Tt{T^'{Q)So{H)11{H)So{H)T^{-Q)So{P) 
+ T^\Q)S^{H)T^{-Q)So{P)nP)So{P)} , 



(247) 



where P = H + Q and Sq is the free relativistic propagator for a nucleon in the nuclear medium in 
eq. (|5l|), which can also be written in the equivalent ways: 

1 



So{K) 



— m? + ie 

eik-kp) 



+ 2TTi9{kF - k)5{K^ - 'm^)e{ko] 



+ 



e{kF - k) 



(248) 



K"^ — m? + ie K"^ — m? — ieko 
The self-energy function S is given by eq. (|40|). 

In order to simplify the calculation of the above polarization propagator, we will simultaneously 
compute the two diagrams contributing to eq. ( ^471) . First we note that eq. (|247] ) can be rewritten 
as: 



n^'g(Q) = n^o^(Q)+n^[(g), 



(249) 



where we introduce n^J^(Q) as the polarization propagator shown in Fig. 34, containing n self-energy 
insertions T,{H) in the hole line and I insertions S(P) in the particle line, i.e., 



U^JiQ) ^ -iTr J ^^r^{Q)[So{H)mrSo{H)r'^{-Q)[So{PMP)fSo{P) , (250) 



where again P = H + Q. From this expression one can derive, as particular cases, the leading- 
order response (n = / = 0, no interaction lines) and the first-order self-energy response (with one 
interaction line, given by eq. ( |249| )). 

Using the nucleon propagator in the medium written in the form in eq. (|24^ ) , the product of 
n + 1 propagators appearing in eq. (pSOf) can be expressed as a derivative of order n according to 



[So{H)i:{H)TSo{H) 

= [{P + m)^{H)r{P + m) 

= [{P + m)^{H)f{P + m) 



eih-kp) 



+ 



eikp-h) 



(ii"2 _ m2 + ie)'^+i (iJ2 _ m2 - ie/io)"+^ 



1 



a=0 



e{h - kp) 



+ 



0{kF - h) 



H"^ — a — m? + ie — a — m? — iehn 
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Figure 34: Diagrammatic definition of the polarization propagator 11^^ for a ph excitation with 
self-energy insertions in the hole and particle lines. Only the forward diagram (a) contributes to 
the electromagnetic responses, while the backward diagram (b) corresponds to a negative value of 
the energy transfer. 
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[{P + m)^{H)f{P + m) 



1 



n\ da"- 



a=0 



+ 2TTie{kF - h)5{H'^ -a- m'^)9{ho 



(251) 



H"^ — a — m? + ie 

where a parameter a, which at the end is going to be zero, has been introduced in the propagator 
denominators. A similar equation holds for the propagation of a particle introducing a second 
parameter (3. The polarization propagator n^J" can then be written as 

dh^Sh 



d" 


s 




da^ 




/ 

/3=0 




1 






— a — m 


2 + ie 




1 





+ 2TTie{kF - h)6{H^ -a- m^)e{ho] 
+ 2me{kF - p)6{P'^ -13- m^)e{po 



_P'^ — (3 — m? + ie 
with P = H + Q, and where we have introduced the functions 



(252) 



li^iiH, P, Q) = li^^iho, h; po, p; go, q) 
= ^Tt {r^(Q)[(/f + m)m)nP + m)T''{-Q)[{P + m)S(P)]'(f + m)} . (253) 

The product of the two brackets inside the integral in eq. (|252| ) gives rise to four terms. The 
first of these contains the product of the two free propagators, namely 

^ ^ (254) 



a — + ie P"^ — (3 — m? + ie 



and yields a genuine vacuum contribution, H^^^^^ {Q), which diverges after integration. Therefore 
we subtract out its contribution, since it pertains to a domain beyond nuclear physics. Performing 
this subtraction of the vacuum propagator we obtain 



n::r(Q)-nif'^(Q)=2^ 









da" 




I 

13=0 •' 



p)5{P^ -(3- m^)6{po) 0{kF - h)5{H^ -a- m'^)e{ho) 



+ 



H"^ — a — m? + ie H"^ — P — + 

+ 2TTie{kF - p)0{kF - h)6{P^ - p- m'^)6{H^ -a- m^)e{Po)e{Ho) 



(255) 



Taking the imaginary part according to eq. (^) we obtain the corresponding hadronic tensor^ 



V 



Im 



vr 



®The extra factor V appears since we are computing the response function of an extended system, see eq. (17.17) 
of ref. il 
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27r 









da" 


a=0 dp'- 


I 

13=0 •' 



^^CiH, P, Q)6{P' -P- n?)5{H^ - a - m^) 



X [e{kF - p)9{po) + 9{kF - h)9{ho) - 2e{kF - p)9{kF - h)e{po)9{ho)] 
Now the factor containing the step functions can be expressed in the form 

[9ikF - p)9ipo) + 9ikF - h)9{ho) - 29{kF - p)9{kF - h)9{po)9{h^)] 
= 9{kF- h)9{K)[l -9{kF- p)9{po)] +9{kF- p)9ipo) [1 - 9{kF - h)9iho)] 

so that the hadronic tensor can be written as a sum of two pieces 



(256) 



-im [nr 

vr 



nl 



n 



(O)Mi'i 
nl J 



nl 



nl 



(257) 



(258) 



where 



= 2ttV 



d" 


d' 




da" 


a=0 dp' 


I 

/3=0 



dhocfh 



li^'i'iH, P, Q)5{P^ -p- rr?)5{H^ - a 



m 



(2^)4 

x9{kF- h)9{ho) [1 -9{kF- p)9{po)] (259) 

corresponds to the hadronic tensor sought for electron scattering (Fig. p4|(a)), whereas the second 
term 

dho(fh 



nl (Q) 

= 27tV 



d" 


d' 




da" 


a=0 dp' 


I 

13=0 



(2^)4 -CiH, P, Q)5{P' -P- m')5{H' - a - 
x9{kF- p)9{po) [1 -9ikF- h)9{ho)] 



(260) 



corresponds to a process with negative energy transfer (Fig. 34(b)); hence it does not contribute 
to the electron scattering response and should be disregarded. 

Finally, the integration with respect to ho in eq. ( |259| ) can be explicitly performed by using the 
(5-functions. One then gets the following expression for the n/-th SE contribution to the hadronic 
tensor 



nl 



V 



d" 


d' 




da" 


dP' 


J 



d^h li:'^iE'^{a),h-E'^iP),p;q) 



(27r)3 



4E'^{a)E'Jp) 



X 5{E^,ia) + qo- E'^{P))9{kF - h)9{p - kF) , (261) 
where p = h + q and we have defined the following energy functions of the parameters a, P 

EM 
E'^iP) 
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yjh? + a + m? 



+ /? + rrfl- 



(262) 
(263) 



Expression ( p6lD is the general equation for which we are searching. It is one of the (n + /)-th order 
contributions to the fuh Hartree-Fock hadronic tensor, which is an infinite sum of ah perturbative 
orders. In the particular case n = / = it gives the well known free (OB) hadronic tensor 

^^^B = W^:; = yj ^^^^^^^^^^'^(^h + go - E^)e{k, - h)9{p - k,) . (264) 

Finally, the hadronic tensor corresponding to one self-energy insertion in the particle or hole lines, 
corresponding to diagrams (g) and (h) in Fig. 6 is given by 

W^E = W^:^ + <r , (265) 
where the n = 1, / = term correspond to the the first-order hole self energy diagram (Fig. 6(h)) 

f d^h /fo''(-E^a),h;£;p,p;g) , i? \a(i. h\Qt u \ 

(266) 



a=0 • 



and for n = 0, / = 1 the first-order particle self-energy diagram (Fig. 6(g)) 

f <fh /oT(i^h,h;ii;;(/3),p;g) ^^^ zr^mw^ Mfl^ t ^ 

(267) 



In the above expressions, after the derivatives with respect to the parameters a and (3 are 
taken, the integral over the hole polar angle cos 9^ can be performed analytically by exploiting the 
5-function. Hence the SE contribution to the hadronic tensor can finally be expressed as a double 
integral. Since the self-energy S involves a triple integral, the contribution to hadronic tensor turns 
out to be a 5-dimensional integral, to be carried out numerically. 

D Renormalized self-energy response using the polarization prop- 
agator 

In Appendix C we computed the first-order self-energy contribution to the polarization propagator 
corresponding to the two diagrams of Fig. 6. The corresponding hadronic tensor splits into the sum 
of the two terms given in eqs. ( |26(: , 267] ) with Fock self-energy insertions in the hole and particle 



lines respectively, and reads 
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da 



d?h Ifo"(K(«)>h;ii;p,p;g) 



r d^h Jo7(j^h,h;ii;;(/?),p; 



6iE^{a) + qo- Ep)9ikF - h)e{p - kp) 



-5{E^ + go - E'{f5))0{kF - h)9{p - kp) , 



(268) 



where p = h + q, and the modified energies for holes and particles have been introduced in 
eqs. (|262|J263 ), with a and (3 being real parameters. Finally the functions I^^ are defined in 
eq. (IH). 

In order to prove the equivalence between the responses computed using the polarization propa- 
gator in eq. ( |268| ) and the result in eq. ( |145|) , obtained using the renormalized current and energies, 
we proceed to perform the derivative with respect to a and /?. For a general function -F(/io) we 
have 



dF{Ei{a)) 



da 



1 



a=0 



2Eh 



dF{ho) 



dho 



(269) 



Hence, interchanging the derivatives and the integral, we get for the hadronic tensor the expression 
r d^h 1 n [l^^(ho,h;Ep,p;q)] 

o(-Eh + qo- Ep)e{kF - h)9{p - kp) 



V 



+ V 
+ V 

- V 



(27r)3 AEi,Ep dho 
d^h 1 



2/io 



ho=E^ 



(27r)3 AEy^E^ 
d^h 1 



lfo''(-Eh,h;£^p,p;g) 



1 d 
2Eh dqo 



KEb +qo- Ep)eikF - h)e{p - kp) 



d 



(27r)3 4Sh^p dpo 
d^h 1 



-^01 (-£^h,h;po,P; 



2po 



5(^h + go - Ep)e{kF - h)e{p - kp) 



Po=Ep 



(27r)3 4Sh^j 



Ioi{Ei„h;Ep,p;q) 



1 d 
2Ep dqo 



6{E^ +qo- Ep)9{kF - h)e{p - kp) 



(270) 



In differentiating the function Iq^ defined in eq. ( 253 ), we first consider the term: 



d 

dpo 



1 

So(p) 
2E„ 



2m 



Er- 



Po=Ep 

(f + m) + 7o(f + m) + (f + m)7o 



Po=Ep 



+ 



2Eh''' ' "~' dpo 
where use has been made of the results 

S(P)(f + m) = So(p)(f + m) 
{P + m){P + m) = 2m(f + m), 



VO=Ep 



(271) 

(272) 
(273) 
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which hold for P^^ on-shell and where So(p) is the eigenvalue of the self-energy for on-shell spinors. 

Next we should compute the derivative of the self-energy 5](P). This function has the general 
structure given in eq. (|4^), and its derivative implies derivatives of the coefficients A, B, and C, 
namely 

aS(P) dA(P) dBiP) dC(P) 



-7 • p + B{P)jo 



Opo dpo dpo dpo 

which must be evaluated for P'^ on-shell. Using again eq. (|273D together with the identity 

we obtain, for P on-shell, 



1 



P L 



dpo 

2 dA{P) 2 dBjP) 2 dC{P) . 
m — p h -fc/p-DolPj 



Po=Ep 



{P + m) 



dpo ^ dpo ^ dpo 
= a{p){P + m) , 

where the definition of the function a(p) in eq. ( |100D has been used. 

Finally, collecting the above results, the derivative in eq. ( |271| ) is found to read 

;f + + m) 



d 

dpo 



1 

So(p) 



Po=Ep 



2Er 



— (f + m) + (f + m) P 



m 



Hence the following expression 



d 
dpo 



I^^{H,P,Q) 

Tv|r'^(Q)(^ + m)r-(-Q) 



^o(p) - "T- ^ «(p) 
Ep 2Ep 2 



Tr <^ r^(Q)(^ + m)r^(-Q)(f + m) 



So(p) joEp -m ^ a(p) 



2K 



yields the derivative of /g\ {H, P, Q) with respect to po and the similar result 



d 
dho 



I^^{H,P,Q) 
Tr|r'^(Q) 



- ho=E^ 

So(h)7o£^h ^ a(h) 



E]^ 2Eh 

+ Tr^r^(Q)(/f + m) 



So(h) 70-Eh - m Q(h) 



Ph 2£'h 



(274) 



(275) 



(276) 



(277) 



(278) 



(279) 
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holds for the derivative of I^q{H, P, Q) with respect to /iq- In addition, with the help of eq. ( 272 ), 
we can write for on-shell momenta 



Tr {r^(Q)(/f + m)T.{H){P + m)^ + m)} 

2mSo(h)Tr {r'^(Q)(^ + m)T^{-Q){P + m)} 



and, as well, 



/oT(^^,^,Q) = 2mSo(p)Tr{r'^(Q)(/f + m)r(-Q)(f + m)} 



(280) 



(281) 



Finally the response functions are found as linear combinations of the diagonal components of the 
hadronic tensor, i.e. W^^. Using the above equations the latter reads 
d^h 1 



V 



(27r)3 AE^E^ 
xTr(r'^(g) 



5{E^ + 



So(h) 70-E'h - m ^ a{\i) 



Eh '^Eh 



Eh 2£'h 
+r'^(Q)(^ + m)r^(-Q) 



Ep)e{kF - h)9{p - kp) 

P + m)T''{-Q){P + m) 
r^(-Q)(f + m) 
f + m) 



So(h) joEh - m ^ a(h) 



So(p) -m ^ a(p) 



Ep 2Ep 



+ V 



+ r^(Q)(^ + m)r''{-Q){P + m) 
d^/i 1 



i:o(p) 70-E'p - m a(p) 
2K 



(27r)3 4^h^p 
m 



^^p 

Tr {r'^(Q)(^ + m)r^(-Q)(/' + m)} 



X So(h) 



So(p) 



m 



d 



Ep I dq^ 



SiEh + qo- Ep)e{kF - h)9{p - kp) 



(282) 



which coincides with the result in eq. (|l4q ), obtained by computing the response functions using 
the renormalized current and energy. 



E On-shell self-energy and field strength renormalization function 



In this Appendix we show in detail how to evaluate the on-shell self-energy in eq. (|82D and the field 



strength renormalization function in eq. (100). They can be expressed in terms of the integrals 
I{P) and L^(P) in eqs. as follows: 



So(p) = 2mB{Ep, p) = -12m^ poLo{po,p) - pL3{po,p) - m^I{po, p) 



Po=Ep 



(283) 



100 



and 



Boip) + 



2dA{po,p) „2^-S(Po,p) 2^C(P0,P) 
m h p 



dpo 



dpo 



dpo 



f2 

mi 



Loipo,p) ^, . , dLo{po,p) p dL^{pQ,Y>) dI{po,p) 



Po 



dpo 



Po dpo 



Po dpo 



PO=Ep 

(284) 



where we have used eqs. (^7[]49|) and the derivatives 



V dpo / p,=Ep 
/ dB{po,p) \ 

V dpo Jp,,=Ep 

fdC{po,p) 



f 



2 r 



Po- 



dLo{po,p) dL3{po,p) 2 9I{Po,p) 



dpo 



dpo 



m 



+Lo(po,p) -PoI{Po,p)] 



mi 



Po- 



Po=Ep 

dLo{po,p) dLs{po,p) 



dpo 



dpo 



m 



\ dpo Jp,=Ep 



-6 



Po- 



dLo{po,p) dL3{po,p) 



+Lo{po,p) 



dpo 
Po 



P 



L3{po,p) 



dpo 



Po=Ep 



m 



dpo 

2dI {po,p) 

dpo 
2dI{po,p) 



Po=Ep 



(285) 
(286) 



dpo 



(287) 



By choosing the z-axis in the direction of p the angular integrals can be performed analytically, 
yielding 

1 



/(^P,P) 

LoiEp,p) 

L3{Ep,p) 
dI{po,p) 



''F k -f{p, k) + 2pk 

(27r)2 Jo ApE], 7(p, k) - 2pk 

1 fkp k , 7(p, A;) + 2pk 
dk— In ■ 



(27r)2 Jo 4p 7(p, k) - 2pk 

1 f A;2 kj{p, k) j{p, k) + 2pk ' 

W)^ Jo \ 2pE^ ~ Sp^Ei, 7(p, k) - 2pk 



dpo 
dLo{po,p) 



1 



dpo 
dL3{po,p) 



dpo 



Po=Ep 
Po=Ep 
Po=Ep 



kp 

dk,^^{Ep-Ei,) 



(27r)2 Jo 2pEk 
1 k 



1 



1 



1 



2p 



"f{p, k) + 2pk 7(p, k) — 2pk 
1 1 



7(j), k) + 2pk 7(p, k) — 2pk 
j{p, k) + 2pk 



(288) 
(289) 
(290) 
(291) 
(292) 



(27r)2 Jo Ap^E], 



j{p, k) — 2pk 



+7(p, k) 

where we have defined the function 



1 



1 



7(p, k) + 2pk 7(p, k) — 2pk 



(293) 



7(p, k) = {Ep - E^f - np^ -k^-ml = 2m^ -ml- 2EpE^ . 



(294) 
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By replacing the above integrals in eqs. ( p83 , ^84[ ) we obtain 



and 



a{p) 



3mf 
2'ir'^m'i 



dk — 



ml 



vr 



p JO 



dk 



^ ^ ^ lip, k) + 2pk 
Apk 7(p, k) — 2pk 

Eli — Ep 



Ek 7^(p, k) — ^p^k"^ 



It is interesting to note that for large p- values the following limits hold 



lim aip) 
lim T,o{p) 

p— >oo 



3m f f^'' ,i,k^ 3m/2 
(XK~ 



TP 1 2^ Ep + kp 

hpkp — mm 



m 



(295) 

(296) 

(297) 
(298) 



where Ep = ykp + m? is the Fermi energy. For kp = 237 MeV/c the on-shell self-energy limit is 
~ 34 MeV. 



F The electromagnetic current operator 



In this Appendix we provide a simple derivation of the non-relativistic reduction of the single- 



nucleon on~shell electromagnetic current operator (see ||66|, [75|, |73[). The single-nucleon electro- 
magnetic current reads 



J^(PV;Ps) = u(p',s' 



Zm 



u{p,s) 



(299) 



where = {E,p) is the four-momentum of the incident nucleon, P'^ = {E',p') the four- 
momentum of the outgoing nucleon and = P'^ — P^ = {u, q) the transferred four-momentum. 
The spin projections for incoming and outgoing nucleons are labeled s and s', respectively. We 



follow the conventions of Bjorken and Drell |31] for the u-spinors. For convenience in the discus- 
sions that follow of the scales in the problem we introduce the dimensionless variables: rj = p/m, 

— For the outgoing nucleon, rj' and 



4m^ 



e = E/m = A = ^, K = 2^ and r 

e' are defined correspondingly. 

For any general operator whose F-matrix form is given by 

Til ri2 

r22 

one has u{p' , s')Fn(p, s) = x^'^Xs, with the current operator F given by 



(300) 



F = ]^^{l+e){l+e') ("fh + Fi2-^ ^ ^ 



T+7 



. 21 



a ■ T] 



(301) 
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An important point in our approach is that we expand only in powers of the bound nucleon 
momentum rj, not in the transferred momentum k, or the transferred energy A. This is a very 
reasonable approximation as the momentum of the initial nucleon is relatively low in most cases, 
since the typical values of rj lie below r]F = kp/ni, where kp is the Fermi momentum {rjp is typically 
about 1/4). However, for those cases corresponding to short-range properties of the nuclear wave 
functions it will be necessary to be very careful with the approximations made. Indeed, for large 
values of r/ a fully-relativistic approach will likely prove necessary. Expanding up to first order in 
powers of rj we get e ~ 1 and e' ~ 1 + 2A. 

Thus, the non-relativistic reductions of the time and space components of the single-nucleon 
electromagnetic current operator can be evaluated in a rather simple form. 

Let us consider first the case of the time component. We have 

J°(P's'; Ps) = n(p', s')j\{p, s) = xl^J^Xs , (302) 

with the current operator = Fi-y^ + iF2a^'^Qu/'2m. Using the general result given by eq. ( ^01| ) 
and expanding up to first order in rj, it is straightforward to get the relation 

^ ^Ge + (Gm -^y>^X'n).a , (303) 

where we have introduced the Sachs form factors Ge = — TF2 and Gm = Fi + F2, and have 
used the relations 

A ~ T + K T] (304) 
~ t{1 + t + 2K-r]) . (305) 
The expression ( |303| ) coincides with the leading-order expressions obtained in previous work ||6^ , 



73|; in those studies a different approach was taken which, while more cumbersome, does yield terms 
of higher order than the ones considered in the present work. It is important to remark again that 
no expansions have been made in terms of the transferred energy and transferred momentum; 
indeed, k, A and r may be arbitrarily large in our approach. 

Let us consider now the case of space components. Thus, we have 

J{P's'; Ps) = n(p', s')Ju{p, s) = x^JXs- (306) 

Introducing the matrix form of the vector component for the single-nucleon electromagnetic current 
operator in the general relation ( |301| ), one can finally write 

J ~ -j^==^^iGM{o- X k) + (^Ge + ^Gai^tj + Gek 
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Gm , s iGE , ^ 

-(k • r])n -\cr x n)n ■ t] 



2(1 + t)' " 2(1 + t) 

- tT{GM - GE/2){a- xrj) + '^^^f'^f^ K x r/)^ • k] , (307) 

2(1 +rj J 

where we have used the relations given by eqs. (|o|,|o|). 



In order to compare with |T^, we write the expression for the transverse component of the 
current, i.e., = J — ^-^n. After some algebra we get the final result 



~ |iG'A/(q- X k) + {ge + i^Gm^ (v - —^^^ 



^Gm ^ _ ^ — ^Gm ^ _ ^ I ^ (308) 



1 + r' ' ' 2(1 +T 
It is straightforward to prove that this expression coincides with the result given by eq. (25) in ||7i 
for an expansion in powers of rj up to first order. 

Therefore, as can be seen from eqs. ( |305 , 30^ ), at linear order in r/ we retain the spin-orbit 



part of the charge and one of the relativistic corrections to the transverse current, the first-order 
convective spin-orbit term. It is also important to remark here that the current operators given by 
eqs. ( |303 , 307l ) satisfy the property of current conservation AJq = k J. Finally, it is also interesting 



to quote the results obtained in the traditional non-relativistic reduction [[7^, |96|, |97|, |98[], where 
it is assumed that k << 1 and A << 1: 

J^nonrel = Ge (309) 



J nonrel = -^G'm['^ X Ct] + G^; 



(310) 



Note that this traditional non-relativistic reduction contains both terms of zeroth and first order 
in r], i.e., the convection current, and is therefore not actually of lowest order in ij. 

We see that the expansion of the current to first order in the variable r] = p/m yields quite 
simple expressions; moreover the various surviving pieces of the relativized current (i.e., charge 
and spin-orbit in the longitudinal and magnetization and convection in the transverse) differ from 
the traditional non-relativistic expressions only by multiplicative (g, a;)-dependent factors such as 
k/\/t or 1/Vl + T, and therefore are easy to implement in already existing non-relativistic models. 
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